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INTRODUCTION 


Background  and  Motivation 

The  need  for  improvements  in  the  state-of-the-art  of  compact,  light- 
weight, low  power,  highly  stable  atomic  frequency  standards  is  exemplified 
by  programs  such  as  that  for  the  NAVSTAR  Global  Positioning  System  (GPS).* 
The  needs  of  the  GPS  program  have  resulted  in  efforts  to  advance  the  state- 
of  -the-art,  especially  for  hydrogen  storage  and  cesium  beam  devices.  These 

needs  include  size,  weight,  and  power  reductions  and  a fairly  stringent 

-13 

long-term  stability  requirement  (1  x 10"  fractional  frequency  stability 
over  periods  of  approximately  one  to  twelve  days  for  phase  II  of  GPS).  Re- 
quirements in  future  programs  of  this  type  can  be  expected  to  be  even  more 
demanding  as  the  overall  state-of-the-art  of  military  technology  moves  for- 
ward. 

Only  three  types  of  frequency  standards  appear  to  have  been  seriously 
considered  for  the  GPS  program:  cesium  beam,  hydrogen,  and  rubidium  gas  cell 
devices.  While  the  cesium  and  hydrogen  devices  have  very  good  long-term 
stability,  they  tend  to  be  bulky,  heavy,  and  to  exhibit  moderate  to  high 
power  consumption.  Rubidium  devices,  on  the  other  hand,  can  be  made  very 
compact,  lightweight,  and  with  low  power  requirements.  (They  also  have  ex- 
cellent short-term  stability.)  As  an  example  of  the  state-of-the-art  in 

2 -3  3 

this  respect,  one  commercial  unit  has  a volume  of  1.1  x 10  m (1.1  liters, 
about  the  size  of  a grapefruit),  weighs  1.3  kg  (2.9  lb),  and  uses  only  13  W 
of  power. 

The  long-term  stability  of  commercial  cesium  devices  is  typically  parts 

13  14  12 

in  10  and  10  over  periods  of  days  to  weeks,  and  parts  in  10  over  much 

longer  periods  of  time.  Rubidium  gas  cell  devices  have  generally  been  con- 
sidered not  to  be  competetive  with  cesium  in  this  area,  but  recent  data  have 

shown  that  this  is  an  over  simplified  point  of  view.  For  example,  the  sta- 

3 4 5 

bility  of  at  least  one  commercial  rubidium  device  has  been  measured  ’ and 

found  to  be  parts  in  1013  and  1014  over  approximately  3 x 104  sec  (8  hours). 

A high  reliability  version  of  this  device,  developed  for  use  in  phase  I GPS 

5 13  14 

satellites,  has  shown  even  better  performance:  parts  in  10  and  10  over 
10^  sec  (30  hours). 

To  summarize,  there  is,  at  the  present  time,  no  single  device  among 
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atomic  frequency  standards  that  simutaneously  provides  state-of-the-art 
capability  in  all  areas  of  interest,  namely,  low  power  consumption,  small 
size  and  weight,  good  short-term  frequency  stability,  and  good  long-term 
frequency  stability.  However,  rubidium  devices  appear  to  come  the  closest 
to  satisfying  all  of  these  requirements:  it  is  much  more  likely  that  the 
long-term  stability  of  rubidium  devices  can  be  significantly  improved,  than 
it  is  that  the  size  and  weight  of  cesium  and  hydrogen  devices  can  be  re- 
duced to  that  of  present  day  state-of-the-art  rubidium  devices. 

Long-Term  Stability  of  Rubidium  Devices 

There  are  at  least  three  possible  factors  that  have  the  potential  of 
degrading  the  long-term  stability  of  a rubidium  gas  cell  device.  They  are: 
the  so-called  "light  shift"  effect,  changes  in  buffer  gas  pressure,  and  the 
so-called  "position-shift  effect".  The  first  two  effects  are  well  known 
and  have  been  recognized  as  possible  problem  areas  since  the  early  days  of 
optical  pumping  research. ^ The  "position-shift  effect"  is  a recent  devel- 
opment^ and  results  from  inhomogeneities  in  the  resonance  cell,  to  be  dis- 
cussed later.  The  object  of  the  work  described  in  this  report  is  to  con- 
struct a rubidium  gas  cell  device  in  which  the  light  shift  is  eliminated  to 
measure  the  resulting  long-term  stability  of  this  device  and  then  to  study 
the  effects  of  the  remaining  parameters  on  the  long-term  stability.  It 
should  also  be  mentioned  that  it  appears  possible  with  this  approach  to  elim- 
inate, or  greatly  reduce,  the  position-shift  effect. 

Pulsed  Optical  Pumping 

The  light  shift  effect  can,  in  principle,  be  eliminated  by  using  the 
technique  of  pulsed  optical  pumping.  In  this  method,  the  pumping  light  and 
microwave  radiation  are  alternately  pulsed  on  and  off  at  a rapid  rate  in 
such  a way  that  the  light  and  microwave  radiations  are  never  simultaneously 
present  in  the  resonance  cell,  i.e.,  there  is  no  temporal  overlap  between 
the  light  pulses  and  the  microwave  pulses.  Since  there  is  no  pumping  light 
present  to  shift  the  energies  of  the  ground  state  hyperfine  levels  of  the  Rb 
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atoms  when  the  microwave  radiation  is  Inducing  the  0-0  hyperfine  transition, 
there  is  no  light  shift. 
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Previous  Work 


It  appears  that  the  idea  of  using  pulsed  optical  pumping  first  orig- 
9 10 

inated  with  Arditi.  In  the  original  experiments,  it  was  demonstrated 
that  there  was  no  light  shift  to  within  the  precision  of  the  measurements 
(one  part  in  1010),  but  the  long-term  stability  was  not  investigated.  To 
the  best  of  our  knowledge,  < he  only  work  on  pulsed  optical  pumping  sub- 
sequent to  that  of  reference  10  appears  in  the  Russian  literature.  The 
Russian  work,  which  is  described  below,  has  not  really  addressed  the  ques- 
tion of  frequency  stability  using  pulsed  optical  pumping;  rather  it  has 
concentrated  almost  exclusively  on  the  physics  of  the  pulsed  optical  pump- 
ing process. 

The  paper  by  Yakobson^  gives  a theoretical  treatment  of  pulsed  pump- 
ing using  a three-level  model.  The  density  matrix  equations  are  solved 
using  a Fourier  series  method  to  obtain  approximate  expressions  for  the  sus- 
ceptibility and  line  shape.  This  is  done  for  three  different  detection 
schemes:  (1)  microwave  detection  of  the  oscillating  magnetization  of  the  Rb 
atoms  (method  used  in  reference  10);  (2)  optical  detection  using  a weak  cw 
auxilliary  light  source  (second  Rb  lamp)in  addition  to  the  (pulsed)  pumping 
light  source;  (3)  optical  detection  using  the  pumping  light  (method  used  by 
The  first  of  four  papers  from  Bazarov's  group,  reference  12,  mostly 
contains  experimental  results  for  an  apparatus  that  uses  detection  scheme 
(2),  above.  Data  are  presented  for  two  different  modes :cw  microwave  field, 
and  pulsed  microwave  field.  In  both  cases  there  is  light-microwave  overlap 
because  the  auxilliary  (probe)  light  source  is  operated  cw.  As  a result 
there  is  a light  shift  present.  Graphs  are  given  of  light  shift  for  values 
of  other  parameters.  FM  line  shapes  are  also  given.  This  work  is  of  little 
interest  in  the  present  context. 

The  original  work  by  Arditi  and  Carver^  showed  that  it  is  possible  to 
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employ  the  time-domain  equivalent  of  separated  oscillating  fields  to  pro- 
duce a "Ramsey  pattern"  for  the  line  shape.  This  is  done  by  using  two  temp- 
orally distinct  microwave  pulses  during  the  dark  time  between  successive 
light  pulses.  The  advantage  of  this  method  is  that  it  allows  the  line  width 
to  be  made  smaller  than  the  natural  line  width  (the  latter  being  character- 
ized by  the  relaxation  time  of  the  Rb  atoms  in  the  gas  cell).  In  these  orig 
inal  experiments,  the  Ramsey  pattern  was  observed  by  microwave  detection  of 
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the  oscillating  magnetization  of  the  Rb  atoms.  The  second  Russian  paper, 
from  Bazarov's  group,  demonstrates  experimentally  that  it  is  also  possible 
to  observe  a Ramsey  pattern  using  optical  detection  with  a weak  auxilliary 
(probe)  light  beam  that  is  pulsed.  In  addition,  this  paper  also  gives  a 
theoretical  expression  for  the  Ramsey  line  shape  based  on  a simplified 
three-level  model. 

15 

The  most  important  Russian  journal  article  is  the  third  paper 
from  Bazarov's  group  which  deals  with  light  shifts  in  pulsed  optical  pump- 
ing. The  pulsed  pumping  scheme  that  is  considered  in  this  paper  is  the 

14 

same  as  that  of  the  previous  paper  from  this  group.  In  this  third  paper, 
the  theory  of  the  previous  paper  is  developed  in  much  greater  detail,  and 
with  many  of  the  previous  restrictions  removed;  experimental  results,  in- 
cluding light  shifts,  are  presented  and  interpreted  using  this  theory.  The 
significance  of  this  paper  is  that  the  light  shifts  under  consideration  are 
occurring  when  there  is  no  overlap  between  the  light  pulses  and  the  applied 
microwave  pulses  (recall  that  the  weak  auxilliary  light  beam  is  also  pulsed, 
as  in  reference  14).  We  believe  that  this  effect  is  very  small  in  our  ex- 
periments, as  is  discussed  later  in  this  report. 

The  last  Russian  paper, ^ also  from  Bazarov's  group,  discusses  (both 
theoretically  and  experimentally)  the  effect  on  the  Ramsey  pattern  of  more 
than  two  microwave  pulses  during  the  dark  period  of  the  pulsed  pumping 
cycle. 

It  is  clear  from  the  Russian  work  that  additional  Russian  work  (not 
described  above)  has  been  done.  Unfortunately  this  latter  work  has  been 
published  in  that  portion  of  the  Russian  literature  that  is  inaccessible  to 
us.  However,  judging  from  the  references  to  this  work  (in  the  papers  that 
are  accessible  to  us),  it  is  of  minor  importance  compared  to  the  work  of 
Bazarov's  group  (references  12,  14-16). 

Present  Work 

The  scope  of  the  work  described  in  this  report  is  shown  in  Figure  1. 

The  design  and  construction  of  the  apparatus  is  discussed  in  the  next  sec- 
tion. Subsequent  sections  of  this  report  deal  with  the  two  other  major 
aspects  of  this  work:  frequency  stability  measurements  (see  Figure  2)  and 
fundamental  device  improvement  (see  Figure  3). 
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Figure  1.  Scope  of  the  work  described  in  this  report. 
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Figure  2.  Flow  chart  for  frequency  stability  me?  <;u  remen  ts ; •=  completed, 
©=  some  work  has  been  done,0=  c anned. 


W.  FUNDAMENTAL 
DEVICE 


Figure  3.  Flow  chart  for  fundamental  device  improvement;  • = completed,  ©=  some  work  has  been  done, 

0=  planned. 


DESIGN  AND  CONSTRUCTION  OF  APPARATUS 


Design  Philosophy 

In  the  original  experiments  of  Arditi  and  Carver,^  the  oscillating 
magnetization  of  the  Rb  atoms  was  detected  using  a sensitive,  microwave 
superheterodyne  receiver.  Besides  adding  considerable  complexity  to  the 
electronic  circuitry,  this  approach  also  requires  use  of  a high  Q microwave 
cavity  for  good  S/N.  This  latter  requirement  has  several  disadvantages  for 
practical  devices:  (1)  small,  low  Q cavities,  which  are  needed  to  reduce 
size  and  weight,  cannot  be  used,  and  (2)  the  device  becomes  susceptible  to 
cavity  pulling.^ 

The  work  of  Bazarov  and  coworkers  has  shown  that  it  is  possible  to  use 
optical  detection  instead  of  microwave  detection.  This  eliminates  the  need 
for  a superheterodyne  receiver,  but  at  the  expense  of  adding  a second  (aux- 
iliary) Rb  lamp  and  its  associated  electronics.  This  is  undesirable  in  a 
compact,  lightweight,  low  power  device  because  the  Rb  lamp  consumes  a sig- 
nificant fraction  of  the  total  power. 

The  design  philosophy  that  has  been  adopted  in  the  present  work  is  that 
it  is  not  profitable  to  conduct  laboratory  experiments  using  principles  and 
methods  that  cannot  subsequently  be  applied  to  the  development  of  a compact, 
light  weight,  low  power  frequency  standard.  For  this  reason,  we  have  chosen 
a substantially  different  approach  to  pulsed  pumping  than  previous  experimen- 
ters. This  has  been  possible,  at  least  in  part,  because  of  our  prior  exper- 
tise in  developing  and  manufacturing  small,  light  weight,  and  low  power  rubi- 
dium frequency  standards  exhibiting  state-of-the-art  performance. 

Our  approach  has  been  to  construct  an  experimental  pulsed  optical  pump- 
ing apparatus  using  as  a starting  point  the  Efratom  Model  FRK  miniature  rubi- 
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dium  frequency  standard  (hereafter  referred  to  as  the  FRK).  The  physics 
package  in  this  unit  (which  includes  the  Rb  lamp  and  associated  electronics) 
is  approximately  cylindrical  in  shape  with  a diameter  of  0.057  m (2.25  in.) 
and  a length  of  0.10  m (4.0  in.).  The  small  size  is  made  possible  mainly  by 
the  substitution  of  a TE^  microwave  cavity  for  the  more  conventional  but 
larger  TEg^  cavity,  and  by  eliminating  the  need  for  a separate  filter  cell. 
The  latter  is  accomplished  by  combining  the  filtering  function  with  the  re- 
sonance function  in  the  resonance  cell.  This  dual  function  cell  is  usually 


referred  to  as  an  "integrated  cell". 

The  pulsed  optical  pumping  apparatus  that  we  have  constructed  uses  an 
unmodified  FRK  physics  package.  The  Rb  lamp  operates  on  the  plasma  dis- 
charge principle  using  a 90  MHz  driving  source.  The  lamp  is  pulsed  on  and 
off  by  pulsing  the  dc  supply  voltage  of  the  90  MHz  lamp  driver. 

In  our  experimental  apparatus  the  Rb  resonance  is  detected  optically 
by  observing  the  absorption  of  the  pumping  light.  This  eliminates  the  need 
for  a sensitive  microwave  receiver  and  a high  Q microwave  cavity.  Our  de- 
tection scheme  also  eliminates  the  need  for  a second  (auxilliary)  Rb  lamp. 
This  is  why  it  has  been  possible  to  use  an  unmodified  FRK,  miniature  physics 
package. 

The  pulsing  scheme  for  these  experiments  is  shown  in  Figure  4.  The 
lamp  is  pulsed  on  for  a time  T^.  After  the  lamp  is  turned  off,  a time  T$ 
elapses  before  the  microwave  pulse  is  pulsed  on.  The  duration  of  the  micro- 
wave  pulse  is  T^.  Following  the  end  of  the  microwave  pulse,  is  a period  of 
time  T^  before  the  beginning  of  the  next  cycle.  All  times  typically  lie 
in  the  range  of  0.1  to  10  msec.  As  can  be  seen  from  the  figure,  the  dark 
time  T(j  = + T^  , and  the  cycle  time  Tc  = T^  + T^.  So  far,  we  have  used 

only  a single  microwave  pulse  between  successive  light  pulses.  We  could 
very  easily  obtain  a Ramsey  pattern  by  modifying  our  pulsing  circuitry  so 
that  our  single  microwave  pulse  of  Figure  4 would  be  replaced  by  two  (or 
more)  microwave  pulses.  We  have  chosen  not  to  do  this  at  the  present  time 
because  of  the  added  complexity  that  it  introduces  when  interpreting  the  ex- 
perimental results;  moreover,  it  is  not  even  certain  that  this  could  produce 
any  improvement  in  device  performance. 

As  can  be  seen  from  Figure  1 (and  Figure  3),  two  parallel  approaches 

have  been  employed.  The  first  involves  making  frequency  stability  measure- 

4 5 

; ments,  usually  over  periods  of  4 x 10  sec  (12  hours)  to  2 x 10  sec  (2H 

days),  and  seeking  to  improve  device  stability  by  means  that  do  not  require 
a fundamental  understanding  of  device  operation.  The  second  involves  a 
fundamental  investigation  of  device  characteristics  for  the  purpose  of  im- 
proving device  operation  as  a frequency  standard.  In  the  first  approach 
(device  operating  as  a frequency  standard)  it  is  necessary  t • frequency  mod- 
ulate the  atomic  resonance  (FM  mode  of  operation).  In  the  second  approach 
it  is  desirable,  in  addition,  to  be  able  to  amplitude  modulate  the  atomic 
resonance  (AM  mode  of  operation)  because  this  makes  certain  device  para- 
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meters  accessible  (e.g.,  atomic  line  shapes  and  relaxation  times)  that 
cannot  be  simply  determined  using  the  FM  mode.  For  this  reason,  two 
different  experimental  configurations  (AM  and  FM)  have  been  used.  Both 
are  described  below. 

AM  Configuration 

Initial  considerations  involving  relaxation  times  and  resonance  line 
widths  indicated  that  the  frequency  of  the  pulsing  cycle  for  optimum  sig- 
nal (AM  or  FM)  would  be  of  the  order  of  several  hundred  hertz.  The  mod- 
ulation frequency,  i.e.,  the  frequency  for  phase  sensitive  detection  ("lock- 
in"  frequency),  appears  as  a small  modulation  of  the  amplitude  of  the  light 
pulses  reaching  the  photodetector.  Since  this  modulation  amounts  to  about 
1(j  to  10"^  of  the  light  pulse  amplitude,  it  is  imperative  to  be  able  to 
separate  the  signal  component  at  the  modulation  frequency  from  the  undesir- 
ed component  at  the  pulsing  frequency.  This  can  be  done  using  analog  fil- 
tering techniques  if  the  modulation  frequency  differs  appreciably  from  the 
pulsing  frequency.  For  this  reason,  the  modulation  frequency  was  chosen  to 
be  10  Hz. 

Figure  5 is  a diagram  of  the  AM  pulsed  optical  pumping  apparatus.  To 
expedite  and  simplify  apparatus  construction,  electronic  circuits  and  com- 
' ponents  which  comprise  the  Model  FRK,  commercial  rubidium  frequency  stan- 

dard were  used  where  possible.  Even  so,  it  was  necessary  to  design  and  con- 
struct some  of  the  blocks  using  discrete  components  and  integrated  circuits. 
These  blocks  included  the  pulsing  electronics  (see  Figure  6 for  details), 
the  10  Hz  signal  conditioner  and  amplifier  (see  Figure  7 for  details),  and 
the  10  Hz  modulation  oscillator.  In  addition,  it  was  necessary  to  build  a 
i 5 MHz,  pulsed  rf  switch  and  buffer  amplifier.  This  device  is  included  in 

the  microwave  synthesizer  block  in  Figures  5 and  8. 

The  AM  operation  is  best  explained  by  reference  to  Figure  5.  Lf  the 
device  is  operated  cw,  the  Rb  lamp  runs  continuously  and  the  6.8  GHz  micro- 
wave  radiation  is  effectively  switched  on  and  off  (square  wave)  at  a 10  Hz 
rate  ("lock-in"  frequency).  When  the  6.8  GHz  is  tuned  to  the  center  of  the 
Rb  hyperfine  resonance  (F=l,  Mp  = 0 — •-  F=2,  Mp  = 0),  the  light  from  the  Rb 
lamp  that  passes  through  the  Rb  resonance  cell,  and  is  incident  on  the 
photodetector,  is  modulated  at  10  Hz.  This  10  Hz  component  is  amplified. 
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Figure  5.  Block  diagram  of  AM  pulsed  optical  pumping  apparatus. 
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Figure  6.  Block  diagram  of  pulsing  electronics. 
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- OUTPUT 


Figure  7.  Block  diagram  of  10  Hz  signal  conditioner  and  amplifier. 


Figure  8.  Block  diagram  of  FM  pulsed  optical  pumping  apparatus  conn- 
ected as  a frequency  standard. 
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synchronously  detected,  filtered,  and  displayed  on  a chart  recorder.  The 
output  from  the  VCXO  is  multiplied  (by  684)  and  synthesized  to  obtain  the 
6.8  GHz.  The  frequency  of  the  6.8  GHz  is  controlled  by  a stable  dc  power 
supply  that  determines  the  frequency  of  the  10  MHz  voltage  controlled  cry- 
stal oscillator  (VCXO). 

Normally,  the  device  is  not  operated  cw;  the  light  and  microwave  rad- 
iation are  pulsed  using  the  scheme  shown  in  Figure  4.  The  AM  operation  is 
then  modified  as  follows.  The  Rb  lamp  is  pulsed  on  and  off  at  the  pulsing 
frequency,  and  is  unaffected  by  the  10  Hz  modulation.  During  half  of  the 
10  Hz  modulation  cycle  the  microwave  radiation  is  turned  off.  During  the 
other  half  of  the  10  Hz  cycle,  the  microwave  radiation  is  pulsed  on  and  off 
at  the  pulsing  frequency  (according  to  Figure  4)  so  that  there  is  no  light- 
microwave  overlap.  Since  the  pulsing  frequency  is  of  order  300  Hz,  there 
are  approximately  15  cycles  of  light  pulsing  that  occur  during  each  half 
of  the  modulation  cycle.  It  should  be  noted  that  the  pulsing  frequency 
and  the  modulation  frequency  are  derived  from  independent  oscillators  so 
that  the  pulsing  is  asynchronous  with  the  10  Hz  modulation. 

In  practice  it  is  undesirable  to  pulse  the  6.8  GHz  microwave  radiation 

completely  on  and  off  because  of  potential  rf  pickup  problems.  Instead, 
the  microwave  radiation  is  effectively  turned  off  by  shifting  its  frequency 
by  ~ 5 MHz.  This  method  works  because  the  5 MHz  shift  is  much  greater  than 

the  width  of  the  atomic  resonance  ( ~ several  kHz  in  the  worst  case).  In 

Figure  5 (and  also  in  Figure  8),  the  6.834...  GHz  is  synthesized  by  multi- 
plying the  VCXO  frequency  by  six  to  obtain  60  MHz.  The  60  MHz  is  then  mul- 
tiplied by  114  and  mixed  with  5.3125  MHz  in  a step  recovery  diode  to  yield 
6.834...  GHz.  In  our  scheme  of  pulsing,  the  5.3125  MHz  is  gated  on  and  off 
using  a pulsed  rf  switch  and  buffer  amplifier. 

The  apparatus  of  Figure  5 has  been  used  to  measure  changes  in  AM  signal 
as  a function  of  device  parameters.  For  most  measurements  of  this  type,  the 
VCXO  is  tuned  to  a frequency  that  corresonds  to  the  Rb  resonance  frequency, 
and  other  parameters  (T^  , T , T^,  etc)  are  then  varied  (one  at  a time).  If 
the  VCXO  frequency  is  varied  while  all  other  parameters  are  held  fixed,  one 
obtains  the  atomic  line  shape  for  those  experimental  conditions. 

To  avoid  possible  confusion,  it  should  be  kept  in  mind  that  AM  opera- 
tion is  always  open  loop;  i.e.,  servo-type  (or  feedback)  operation  is  not 
possible  using  the  AM  configuration  shown  in  Figure  5. 
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FM  Configuration 


The  FM  configuration  shown  in  Figure  8 allows  the  pulsed  pumping 
apparatus  to  be  operated  as  a frequency  standard.  The  FM  method,  as 
used  here,  consists  of  sinusoidally  sweeping  the  frequency  of  the  6.8  GHz 
radiation  back  and  forth  across  the  Rb  resonance  at  a 10  Hz  rate.  This 
produces  a 10  Hz  modulation  of  the  light  pulses  reaching  the  photodetector 
which,  when  amplified  and  synchronously  detected,  gives  a dc  output  vol- 
tage as  shown  schematically  in  Figure  9.  The  physics  package,  and  this 
portion  of  the  electronics,  therefore  function  as  a frequency  discrimin- 
tor  that  gives  a negative  (positive)  output  voltage  when  the  frequency  of 
the  6.8  GHz  microwave  radiation  lies  above  (below)  the  Rb  resonant  fre- 
quency. If  this  error  signal  is  fed  back,  and  used  to  control  the  VCXO 
frequency,  operation  as  a frequency  standard  results. 

The  FM  configuration  of  Figure  8 differs  from  the  AM  configuration  of 
Figure  5 in  the  following  respects: 

The  10  Hz  modulation  oscillator  is  no  longer  used  to  pulse  the  micro- 
wave  radiation  on  and  off  (effectively).  Instead,  it  is  used  to  phase  mod- 
ulate the  10  MHz  output  of  the  VCXO.  This  produces  frequency  modulation  of 
the  6.8  GHz.  The  advantage  of  phase  modulation  over  frequency  modulation 
is  that  the  frequency  of  the  VCXO  remains  spectrally  pure  (if  frequency 
modulation  of  the  10  MHz  at  a 10  Hz  rate  were  used,  the  VCXO  frequency  spec- 
trum would  contain  sidebands  spaced  at  10  Hz  intervals).  Phase  modulation 
is  accomplished  using  a PIN  diode  phase  modulator.  This  device  has  the  ad- 
vantage of  linearity  for  phase  deviations  that  would  produce  nonlinearity 
if  a varactor  diode  phase  modulator  were  used  instead.  The  PIN  diode  acts 
as  an  electrically  - controllable  rf  resistor  and  is  used  in  a standard  RC 
phase  shifter  configuration.  The  optimum  value  of  the  peak  phase  deviation 
at  10  MHz  has  been  determined  experimentally  to  be  ~3°.  This  corresponds 
to  a peak  frequency  deviation  at  6.8  GHz  that  is  somewhat  greater  than  the 
Rb  resonance  line  width. 

In  addition  to  the  phase  modulator,  a 10  Hz  variable  phase  shifter  was 
designed  and  built  so  that  the  10  Hz  signal  could  be  phase  shifted  relative 
to  the  10  Hz  reference.  This  is  required  for  the  phase  detection  scheme 
used.  The  synchronous  switch  and  integrator  were  taken  from  an  FRK  and 
modified  for  use  here. 
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9.  Schematic  diagram  of  dc  output  of  phase  detector  as  a function  of  microwave  frequency 
configuration  (Figure  8,  with  loop  open),  f is  the  Rb  resonance  frequency,  6. 834... GHz 


FREQUENCY  STABILITY  MEASUREMENTS 


As  mentioned  previously,  two  parallel  approaches  have  been  used 
in  this  work.  The  first,  which  is  discussed  here,  consists  of  making 
frequency  stability  measurements  and  seeking  to  improve  device  stability 
by  means  that  do  not  require  a fundamental  understanding  of  device 
physics. 


Methods  of  Measurement 

Two  methods  of  measuring  frequency  stability  have  been  used  in  this 
work.  In  the  first,  the  10  MHz  output  of  the  VCXO  is  divided  by  two  and 
fed  to  the  input  of  a Tracor  Model  527A  frequency  difference  meter  (FDM). 
The  output  of  the  Tracor  is  a dc  voltage  that  is  proportional  to  the 
difference  between  the  VCXO  frequency  and  that  of  a 5 MHz  reference.  The 
reference  ("lab  frequency  standard"  in  Figures  5 and  8)  is  essentially  an 
Efratom  Model  FRK-H  rubidium  frequency  standard  that  has  been  specially 
selected  and  modified  to  provide  superior  performance  and  temperature  sta- 
bility. The  dc  output  of  the  FDM  is  passed  through  an  RC  smoothing  net- 
work and  then  displayed  on  a chart  recorder. 

The  second  method  uses  a custom-built  TTL  phase  comparator.  This 
device  gives  a dc  output  that  is  proportional  to  the  phase  difference 
between  the  VCXO  frequency  and  the  reference  ("lab  standard")  frequency, 
modulo  2tt. 

The  FDM  is  more  convenient  to  use,  but  it  is  difficult  to  relate 

its  output  indication  to  the  commonly  accepted  measure  of  frequency  sta- 

19 

bility,  the  Allan  variance.  Also,  the  RC  time  constant  is  limited  to 
about  100  sec.  On  the  other  hand,  the  Allan  variance  can  be  determined 
for  arbitrarily  long  averaging  times  from  the  phase  comparator  data. 

Problems  and  Solutions 

Figure  2 shows  the  problems  that  have  been  encountered  in  obtaining 
good  frequency  stability  with  the  pulsed  pumping  apparatus,  and  the 
progress  made  in  finding  solutions. 


Erratic  Frequency  Jumps 

Considerable  time  was  spent  trying  to  eliminate  spontaneous  fre- 
quency jumps.  Typically,  these  jumps  were  several  parts  in  lO1^  and 
would  appear  as  a step  function  in  frequency.  This  problem  was  par- 
ticularly difficult  to  solve  because  it  was  due  to  several  causes. 

Three  methods  of  attack  were  utilized  in  solving  this  problem.  The 
first  was  the  improvement  of  the  rf  pulsing  scheme.  This  was  carried 
out  by  designing  and  building  a 5.3  MHz  buffer  amplifier  and  pulsed  rf 
switch.  This  circuit  pulses  the  5.3  MHz  on  and  off  to  provide  milli- 
second pulses  of  5.3  MHz  (with  rise  and  fall  times  of  < 10  usee)  that 
are  fed  to  the  step  recovery  diode.  The  first  stage  is  an  emitter 
follower  buffer  that  isolates  the  5.3  MHz  source  from  the  rf  switch.  A 
PIN  diode  (HP  3080)  rf  switch  follows  this  stage.  The  rf  pulses  from 
the  PIN  diode  network  are  amplified  and  cleaned  up  using  a tuned  FET 
amplifier,  and  then  passed  into  the  final  stage  which  is  an  emitter 
follower  driver.  The  previous  scheme  for  pulsing  the  5.3  MHz  was  much 
simpler  but  the  pulses  were  not  clean,  and  the  5.3  MHz  source  was  mod- 
ulated at  the  pulsing  frequency.  The  present  scheme  eliminates  these 
problems. 

The  second  method  of  attack  was  much  simpler,  but  at  least  as 
effective.  The  pots  used  to  adjust  the  pulsing  parameters  (T«  , T 
and  Tc;  see  Figure  6)  were  replaced  by  fixed  metal  film  resistors.  This 
eliminated  small  jumps  in  the  pulse  parameters  due  to  small  changes  in 
pot  resistances.  Related  to  this  was  the  third  method  of  attack,  which 
involved  selecting  values  of  the  pulse  parameters  to  reduce  the  sensi- 
tivity of  the  output  frequency  to  small  changes  in  these  parameters. 

Since  these  measures  have  been  instituted,  frequency  jumps  have 
ceased  to  be  a problem. 

Dc  Offsets 

Dc  offsets  in  the  electronics,  particularly  when  they  occur  in  the 
signal  channel  between  the  photocell  and  the  VCXO  (see  Figure  8),  pro- 
duce corresponding  offsets  in  the  output  frequency.  As  an  example  of 
this,  consider  what  happens  when  a dc  offset  is  applied  to  the  input  of 
the  integrator  in  Figure  8.  This  dc  voltage  causes  the  output  voltage 
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of  the  integrator  to  begin  changing.  This,  in  turn,  changes  the  VCXO 
frequency  which  detunes  the  6.8  GHz  from  the  Rb  resonance  frequency. 

This  detuning  generates  a dc  voltage  at  the  output  of  the  phase  detect- 
or (in  accordance  with  Figure  9)  which  tends  to  oppose  the  dc  offset 
that  was  originally  applied.  The  servo  loop  therefore  seeks  to  adjust 
itself  in  such  a way  as  to  produce  a null  condition  at  the  integrator 
input  (net  dc  input  offset  = 0).  Obviously,  the  only  way  that  it  can  do 
this  is  to  detune  the  VCXO  (and  therefore  the  6.8  GHz)  sufficiently  to 
produce  an  error  voltage  at  the  input  of  the  integrator  that  just  cancels 
the  original  dc  offset  (the  dc  loop  gain  for  our  loop  is  ~10^  so  that 
finite  gain  effects  are  unimportant  here).  A detailed  analysis  shows 
that  the  frequency  offset  is  proportional  to  the  dc  offset,  to  the  gain 
of  the  preceding  electronics,  and  to  the  gain  of  the  physics  package 
(all  suitably  defined).  The  output  frequency  therefore  becomes  dependent 
on  the  gain  of  the  electronics  and  the  physics  package  which  is  clearly 
very  undesirable  since  it  cannot  be  assumed  that  these  gains  remain 
constant,  especially  over  long  periods  of  time. 

Another  difficulty  occurs  when  the  dc  offset  is  due,  for  example, 
to  rf  pickup  (mostly  60  MHz  from  the  microwave  synthesizer  and  the  physics 
package,  wherein  is  located  the  step  recovery  diode).  One  way  this  can 
occur  is  if  there  is  partial  rectification  of  the  rf  due  to  nonlinearity 
in  the  signal  electronics.  This  type  of  dc  offset  is  particularly  trouble 
some  because  it  is  difficult  to  eliminate,  and  it  is  sensitive  to  the 
relative  positions  of  circuits,  wires  and  ground  locations,  all  of  which 
are  subject  to  change.  This  effect  was  eliminated  in  our  apparatus  by 
improving  the  rf  shielding,  grounding,  and  decoupling.  A significant 
effort  was  necessary  to  accomplish  this. 

A somewhat  related  problem  is  the  presence  of  dc  offsets  due  to 
ground  loops,  and  inadequate  power  supply  decoupling.  With  care,  these 
problems  can  be  essentially  eliminated. 

All  of  the  offsets  discussed  so  far  can  be  easily  detected  by 

varying  the  gain  of  the  electronics.  On  the  basis  of  this  type  of 

measurement,  it  is  estimated  that  dc  offsets  of  the  type  discussed  can 

be  reduced  to  the  extent  that  a 10%  change  in  gain  produces  a frequency 
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shift  of  less  than  1 x 10  . (The  actual  shift  may  be  much  lower  than 
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this;  1 x 10  is  an  upper  limit  imposed  by  the  resolution  of  our 
measurements. ) 

Another  possible  source  of  dc  offset  is  second-harmonic  contam- 
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ination  of  the  10  Hz  modulation.  ’ If  this  offset  is  constant  there 
is  no  detrimental  effect,  but  this  cannot  be  taken  for  granted.  When 
this  offset  is  present,  changes  in  the  second-harmonic  content  of  the 
modulation  can  produce  output  frequency  shifts.  Even  if  the  amount  of 
second  harmonic  does  not  change,  output  frequency  shifts  can  still  re- 
sult from  changes  in  the  relative  phases  of  signal  and  reference  at  the 
20 

phase  detector  due  to  the  presence  of  second  harmonic  in  the  modulation. 

Second-harmonic  contamination  was  investigated  in  considerable 

detail.  Since  the  problem  is  difficult  to  treat  theoretical ly  for  large 

angle  modulation  (~  10n  at  6.8  GHz),  a circuit  was  designed  to  inject 

known  amounts  of  second  harmonic  into  the  10  Hz  modulation,  so  that  the 

corresponding  shifts  in  the  10  MHz  output  frequency  could  be  observed. 

Using  this  method,  it  was  found  that  the  shift  in  the  output  frequency 

is  proportional  to  the  amplitude  of  the  in-phase  component  (i.e., 

cos  2co_t  for  cos  w t modulation,  where  m/  2tt  = 10  Hz)  of  second  harmonic, 
m m m 

It  was  also  verified  (as  expected  from  theoretical  considerations)  that 

the  frequency  offset  is  independent  of  the  gain  of  the  signal  electronics, 

which  makes  it  difficult  to  detect  its  presence.  For  cw  operation  of 

1 the  frequency  standard,  an  in-phase  second-harmonic  contamination  of  -23 

.0 

dB  produced  a frequency  offset  of  1 x 10  . (The  result  for  pulsed  oper- 

ation is  not  expected  to  be  significantly  different.) 

Measurements  were  also  made  of  the  second-harmonic  content  of 
the  10  Hz  analog  signal  used  to  drive  the  10  MHz  PIN  diode  phase  modu- 
lator. To  do  this,  a simple,  low-frequency  Fourier  analyzer  was  con- 
i' structed.  An  upper  limit  of  -66  dB  was  established  for  the  amount  of 

second  harmonic.  If  this  amount  of  in-phase  second  harmonic  contamina- 
tion were  present,  the  resultant  frequency  offset  would  be  approximately 
-12 

7 x 10  . It  is  likely  that  the  actual  offset  is  considerably  less 

than  this. 

Even  if  the  10  Hz  driving  signal  contains  no  20  Hz  component,  it 
is  still  possible  to  generate  a 20  Hz  component  in  the  modulation  cir- 
cuitry (due  to  possible  nonlinearity  of  the  modulator).  Unfortunately 
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there  appears  to  be  no  simple,  direct  way  to  investigate  this  possibil- 
ity. An  indirect  method  consists  of  observing  the  amount  of  out-of- 
phase first  harmonic  in  the  signal  from  the  physics  package  (which  is 
due  to  the  out-of-phase  component  of  second  harmonic  in  the  modulation) 
when  the  servo  loop  is  closed  and  the  VCXO  is  locked  to  the  Rb  resonance. 
We  have  done  this  for  our  apparatus  and  have  found  negligible  out-of- 
phase first  harmonic.  This  strongly  suggests  (but  does  not  prove)  that 
there  is  negligible  second  harmonic  contamination  of  the  modulation  in 
our  apparatus. 

It  should  be  mentioned  that  several  modifications  of  the  modulation 
circuitry  were  made  prior  to  the  above  measurements.  The  10  Hz  analog 
PIN  diode  driver  was  redesigned  to  provide  more  distortion-free  operation, 
and  the  10  Hz  amplitude  and  dc  bias  controls  were  made  orthogonal . Most 
importantly,  the  PIN  diode  modulator  itself  was  redesigned  to  use  five 
PIN  diodes  (HP  3080)  rather  than  a single  diode.  This  allows  the  use  of 
higher  drive  levels  and  reduces  the  susceptibility  to  noise  and  spurious 
signals. 

To  summarize,  modifications  of  the  modulation  circuitry  have  been 
made  that  have  resulted  in  improved  operation.  Measurements  made  on  this 
improved  circuitry,  and  on  the  apparatus  as  a whole,  have  indicated  that 
dc  offsets  and  frequency  offsets  due  to  second-harmonic  contamination  of 
i the  10  Hz  moaulation  are  small,  if  they  exist  at  all,  and  that  degradation 

of  frequency  stability  due  to  these  effects  is  probably  negligible. 

Short-Term  Stability  (t  = 1 to  100  sec) 

1 

The  signal-to-noise  ratio  of  the  apparatus  determines  the  short-term 
frequency  stability  for  averaging  times  greater  than  the  attack  time  of 
the  servo  loop  (~  0.2  sec).  It  has  been  found  that  the  signal  in  the 
pulsed  mode  is  quite  good:  generally  about  20%  to  30%  of  that  obtained  in 
the  cw  mode.  The  noise,  however,  is  significantly  greater  in  the  pulsed 
mode.  This  is  not  surprizing  since  the  10  Hz  signal  component  is  of  the 

_3 

order  of  10  of  the  amplitude  of  the  pulses.  If  the  light  pulses  exhibit 
spurious  (noise)  amplitude  changes  from  pulse  to  pulse,  or  if  there  is 
light-pulse-width  jitter,  then  the  10  Hz  Fourier  components  will  produce 
frequency  noise.  Clearly,  changes  in  amplitude  and  pulse  duration  as 
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small  as  10  have  the  potential  of  producing  10  Hz  noise. 

It  has  been  found  possible  to  reduce  the  10  Hz  noise  by  a factor  of 
~ 7 by  using  a "light  compensation"  scheme  which  employs  two  photodetec- 
tors. Besides  the  usual  photodetector,  a second  one  is  placed  in  front 
of  the  Rb  resonance  cell,  close  to  the  Rb  lamp  (see  Figure  8).  The  two 
detectors  are  connected  so  that  their  output  pulses  tend  to  cancel. 

Since  only  the  usual  photodetector  receives  the  10  Hz  signal  component, 
the  10  Hz  signal  is  unchanged.  This  scheme  greatly  reduces  the  back- 
ground signal.  Using  this  method,  one  typically  observes  a peak-peak 
fractional  frequency  fluctuation  of  1 x 10’*1  with  a 100  second  time  con- 
stant on  the  output  of  the  FDM.  This  is  very  nearly  equal  to  the  frac- 
tional frequency  fluctuation  obtained  in  the  cw  mode  for  a 10  Hz  modula- 
tion frequency.  Additional  measurements  for  1 and  10  second  time  con- 
stants show  that  the  fractional  frequency  stability  is  improving  in  pro- 

portion  to  t 2 where  t is  the  time  constant.  This  shows  that  we  are  deal- 
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ing  with  white  frequency  modulation  noise,  as  is  usually  the  case  for 
passive  rubidium  devices. 

As  noted  previously,  the  FDM  readings  are  not  easily  related  (ana- 
lytically) to  the  usual  measure  of  frequency  stability,  which  is  the 
square  root  of  the  Allan  variance  (denoted  by  a y (t)  where  x is  the  aver- 
aging time).  However,  experience  gained  at  Efratom  with  our  commercial 
units  has  shown  that  o Y (T)  for  x = t can  be  estimated  quite  reliably  from 
FDM  readings  by  dividing  the  FDM  peak-peak  fractional  frequency  fluctua- 
tions by  six  (dividing  by  five  gives  a more  conservative  result).  Since 
our  FDM  readings  for  the  pulsed  pumping  apparatus  have  the  expected  t 2 
dependence,  this  procedure  should  also  be  valid  here.  We  therefore  con- 
clude that  for  our  pulsed  pumping  frequency  standard,  the  short-term  sta- 
Dility  is  given  by 

oY(x)  = 2 x ID'11  t-*5,  1 < x s 100  sec.  (1) 

-12 

In  particular,  for  x = 100  sec,  oy  = 2 x 10  . This  value  is  listed  in 

Table  1 along  with  values  of  Oy  (100  sec)  for  other  frequency  standards. 
Note  that  the  short-term  stability  of  the  pulsed  pumping  device  is  better 
than  that  of  a commercial  cesium  beam,  and  comparable  to  commercial  rubi- 
diums. 
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TABLE  1.  SHORT-TERM  FREQUENCY  STABILITY 
OF  SEVERAL  FREQUENCY  STANDARDS  FOR  A 
100  SEC  AVERAGING  TIME 


Frequency  standard  Fractional 

frequency  stability3 


Pulsed  optical  pumping  (present  work)b 
Hewlett-Packard  5062C  cesium  beam0 
Efratom  FRK-L  rubidium0 
FRK-H  rubidium0 

Oscil loquartz  B-5400  high  stability  quartz0 
Efratom  FRK-H  rubidium,  Serial  No.  2610d 


2 x 10 
7 x 10 

3 x 10 
1 x 10 
1 x 10 
3 x 10 


-12 

-12 

-12 

-12 

-12 

-13 


aSquare  root  of  Allan  variance,  av(x),  for  x = 100  sec. 

L.  I 

•See  text  for  explanation. 

c 

Manufacturer's  specification 

dData  taken  at  Rockwell  International,  Autonetics  Division, 

Anaheim,  California,  17  Nov.  77.  Oy(x  = 100  sec)  = (2.95  ± 0.26)  x 
10  ^ for  FRK-H,  S/N  2610  using  HP  5065  reference,  rounded  to  3 x 10"^. 
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Tests  have  been  conducted  to  determine  the  sources  of  noise  in  the 
pulsed  pumping  apparatus.  Open-loop  noise  measurements  have  been  made, 
f and  the  open-loop  gains  of  the  signal  electronics  and  physics  package  have 

been  determined  as  well.  From  these  data  it  is  a simple  matter  to  calcu- 
late the  contributions  of  the  various  noise  sources  to  the  observed  Qy (t  ) . 

> It  is  found  that  the  lamp  is  the  major  contributor.  However,  if  there 

were  no  lamp  noise, Oy  would  improve  by  only  a factor  of  two.  As  a check 

on  this  method,  the  total  nv  was  calculated  from  the  total  open-loop  noise 

’ 12 

and  found  to  be  (ay)caic  = 1*6  x 10  . The  value  of  ®y  measured  with 

the  servo  loop  closed  is  (<3y)  mgas  = 2 x 10-1^.  (Both  are  for  x = 100 

sec.)  The  agreement  lies  within  the  experimental  uncertainties. 

The  implication  of  the  above  result  is  that  any  effort  to  increase 

S/N  ratio  by  more  than  a factor  of  two,  must  consider  both  lamp  noise  and 

amplifier  noise.  In  fact,  it  appears  possible  to  improve  the  S/N  by  any 

of  several  methods.  One  method  is  to  improve  the  signal.  Under  certain 
87 

conditions  Rb  lamps  can  give  larger  signals  then  we  presently  obtain 
with  the  present  lamp  (that  contains  both  isotopes  of  Rb).  Optimization 
of  lamp  and  cell  temperatures  (simultaneously  with  the  pulsing  parameters) 
could  also  yield  greater  signal.  At  the  present  time,  systematic  studies 
of  this  nature  have  not  been  carried  out  so  that  additional  improvements 
are  probably  possible. 

The  other  method  of  improving  the  S/N  is  to  decrease  the  noise.  It 

i 

may  be  possible  to  reduce  the  lamp  noise  by  proper  choice  of  lamp  operating 
conditions  (tests  have  shown  that  the  noise  is  not  due  to  the  pulsing  elec- 
tronics). Very  little  effort  has  been  expended  in  this  area.  It  also  ap- 
pears to  be  possible  to  reduce  amplifier  noise  at  10  Hz  by  using  a differ- 

Ient  op  amp  for  the  input  stage  of  the  photocell  amplifier.  We  base  this 
statement  on  some  independent  IR&D  work  that  we  have  done  using  our  com- 
mercial units. 

In  summary,  we  note  that  the  presently  measured  short-term  frequency 
stability  of  our  pulsed  pumping  frequency  standard  is  quite  acceptable. 

For  averaging  times  of  1 sec  to  100  sec,  it  exceeds  that  of  the  Hewlett- 
Packard  5061A  cesium  beam  frequency  standard,  and  is  comparable  to  that  of 
commercial  rubidiums  (see  Table  1 for  data  for  x = 100  sec).  Nevertheless, 
it  appears  possible  to  improve  the  short-term  stability  further,  perhaps 
by  a factor  of  two  or  more,  if  additional  effort  is  expended  in  this  area. 
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Long-Term  Stability  (t  = 24  hours) 


Measurements  of  the  long-term  stability  of  the  pulsed  pumping 
apparatus  have  been  greatly  hampered  by  frequency  drifts  caused  by 
changes  in  ambient  temperature.  There  are  two  basic  difficulties.  The 
first  is  that,  at  the  present  time,  the  apparatus  is  more  sensitive  to 
temperature  changes  than  the  commercial  unit.  This  problem  is  consid- 
ered further,  below.  The  second  problem  is  that  in  the  room  where  the 
apparatus  was  temporarily  set  up,  there  is  a 24  hour  temperature  change 
of  about  7 °C.  At  the  time  of  this  writing,  the  apparatus  has  been  moved 
to  a more  benign  location  (where  the  temperature  variation  is  only  1 to 
2 °C)  and  measurements  are  in  progress  there  using  the  house  frequency 
standard  as  a reference. 

In  spite  of  the  difficulties  encountered  with  temperature-induced 
frequency  changes,  it  has  still  been  possible  to  obtain  some  preliminary 
results  on  long-term  stability  and  drift.  The  method  used  is  to  take  phase 
comparison  data  over  a period  of  a couple  of  days.  The  data  show  signifi- 
cant frequency  fluctuations  that  are  due  to  temperature,  but  these  varia- 
tions are  quasi  periodic  with  a period  of  24  hours,  and  are  correlated  with 
the  24  hour  ambient  temperature  cycle.  If  the  temperature  cycle  were  exac- 
tly periodic  and  if  the  apparatus  were  in  a true  steady-state  condition  in 
response  to  ambient  temperature  changes,  then  any  changes  in  the  average 
frequency  from  one  24  hour  period  to  the  next  would  be  due  to  long-term  fre- 
quency instability  and  drift.  If  the  above  conditions  are  not  satisfied, 
then  additional  frequency  changes  due  to  temperature  effects  would  be  super- 
imposed on  the  long-term  stability  and  drift. 
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Phase  comparison  data  were  taken  over  a weekend,  from  6:30  pm  on  Fri- 
day until  8:00  am  on  Monday.  The  maximum  frequency  (t  = 20  min)  variation 
due  to  temperature  changes  during  this  period  was  approximately  4 x 10"11. 

The  phase  comparison  data  were  analyzed  at  24  hour  intervals  at  times  of  6 am 
and  11  pm.  (These  times  of  day  were  chosen  because  the  frequency  then  was 
changing  slowly  due  to  temperature-)  Table  2 summarizes  the  results  obtained 
for  the  average  frequencies  over  these  24  hour  periods.  It  now  remains  to 
consider  the  effect  on  these  results  of  deviations  from  the  ideal  temperature 
cycle  assumed  above. 

Measurements  have  been  made  of  the  weekend  temperature  variation  in  the 
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room  where  the  stability  data  were  taken.  These  measurements  showed 
a -1°C  temperature  change  from  11  pm  Friday  to  11  pm  Saturday.  All  24 
hour  changes  occuring  from  Saturday  onward  were  also  negative,  with  mag- 
nitude s 0.5  °C.  Since  the  overall  temperature  coefficient  of  the  ap- 
paratus  during  the  stability  run^J  was  about  5 x 10"  /°C,  we  can  there- 

fore conclude  that  the  effect  of  the  deviation  from  a perfect  24  hour 

temperature  cycle  is  to  introduce  a small  negative  frequency  drift  of  < 5 
-12 

x 10  which  is  greatest  early  in  the  run.  Thus,  the  frequency  changes 

due  to  temperature  effects  are  expected  to  be  approximately  the  same  sign 

and  magnitude  as  the  changes  in  average  frequency  that  are  recorded  in 

Table  2.  It  is  therefore  concluded  that  changes  in  the  average  frequency 

12 

(t  = 24  hours)  over  successive  24  hour  periods  are  of  order  parts  in  10 
Since  the  reference  used  (essentially  an  FRK-H)  cannot  be  verified  (by 

-12 

measurements  traceable  to  NBS)  to  be  stable  to  better  than  about  6 x 10 
for  t = 24  hours,  we  have  to  conclude,  on  the  basis  of  these  data,  that 
the  long-term  stability  and  drift  (t  = 24  hours)  for  the  pulsed  optical 
pumping  apparatus  are  no  worse  than  this.  In  fact,  the  results  obtained 
here  do  not  exclude  the  possibility  that  the  stability  of  the  pulsed  de- 
vice is  better  than  this,  and  even  better  than  that  of  the  commercial  fre- 
quency standard.  Please  see  NOTE  ADDED  IN  PROOF,  p.68. 


TABLE  2.  LONG-TERM  FREQUENCY  STABILITY  OVER  24  HOUR  PERIODS 


Twenty-four  hour  period 


Average 
frequency 
(t  = 24  hours) 


Change  in 
average 
frequency 


11  pm  Fri.  to  11  pm  Sat. 
11  pm  Sat.  to  11  pm  Sun. 
6 am  Sat.  to  6 am  Sun. 
6 am  Sun.  to  6 am  Mon. 


-1.7  x 10" 
-2.2  x 10" 
-2.1  x 10" 
-2.3  x 10" 


-5  x 10" 


-2  x 10" 


Additional  long-term  stability  measurements  need  to  be  made  to 
characterize  adequately  the  stability  over  periods  of  24  hours  and 
longer.  The  sensitivity  of  the  output  frequency  to  changes  in  ambient 
temperature  is  an  impediment  to  such  measurements,  and  some  effort  should 
be  devoted  to  reducing  this  sensitivity.  Possible  ways  of  doing  this  are 
discussed  in  the  next  section. 

Measurements  of  long-term  stability  to  better  than  a few  parts  in 

12 

10  will  require  that  a cesium  beam  be  used  as  the  reference:  our  pre- 
sent in-house  frequency  standard  is  a three-member  ensemble.  The  first 
two  members  are  especially  selected  and  modified  Efratom  FRK-H  rubidium 
frequency  standards.  The  third  member  is  the  channel  7 (Los  Angeles)  ABC 
TV  color  subcarrier  which  is  locked  to  a cesium  standard  at  the  network 

originating  studio.  Over  a 24  hour  period,  the  measurement  uncertainty 

24  -12 

for  the  channel  7 color  subcarrier  frequency  is  6 x 10  . For  the  en- 

-12 

semble  as  a whole,  the  long-term  stability  is  64  x 10  for  24  hours 
(8.6  x 104  sec),  and  <2  x 10”12  for  one  week  (6.0  x 10^  sec). 

Temperature  Sensitivity  of  Output  Frequency 

There  are  at  least  three  different  ways  to  reduce  the  sensitivity 

of  the  output  frequency  to  changes  in  ambient  temperature:  (1)  reduce  the 

/ 

temperature  sensitivity  of  the  electronics;  (2)  reduce  the  resonance  cell 
temperature  coefficient;  (3)  use  electronic  temperature  compensation  that 
senses  temperature  changes  and  makes  a correction  to  the  output  frequency. 

The  temperature  sensitivity  of  the  electronic  circuitry  has  been  in- 
vestigated by  cooling  circuits  and  components  with  spray  cooler  (i.e., 
cooling  with  freon  from  a pressurized  can)  while  the  apparatus  is  operating 
as  a frequency  standard,  and  observing  any  associated  change  in  the  output 
frequency.  While  a comprehensive  study  has  not  been  carried  out  using  this 
method,  it  is  still  possible  to  draw  certain  conclusions.  The  5 MHz  buffer 
amplifier  and  rf  switch,  and  the  pulsing  electronics  are  relatively  insen- 
sitive to  temperature  changes  (they  were  designed  to  accomplish  this  end). 
On  the  other  hand,  the  PIN  diode  phase  modulator  is  sensitive  to  tempera- 
ture, especially  the  five  PIN  diodes  in  this  circuit.  Steps  should  there- 
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fore  be  taken  to  either  redesign  this  circuit  for  reduced  temperature 
sensitivity,  or  to  place  it  in  a temperature-controlled  oven  (a  possible 
experimental  expediency).  In  addition,  a systematic  spray  cooling  study 
of  all  other  circuitry  should  be  undertaken,  and  corrective  action  taken 
where  required. 

The  dependence  of  the  output  frequency  on  the  temperature  of  the  Rb 
resonance  cell  has  been  measured  over  a 16  °C  temperature  range. ^ The 
temperature  coefficient  is  approximately  +6.5  x 10"^/°C.  This  is  about 
50%  greater  than  for  our  commercial  unit.  This  temperature  coefficient 
could  be  reduced  by  up  to  10  times  by  using  a different  combination  of 
buffer  gases  in  the  Rb  resonance  cell. 

Electronic  temperature  compensation  (a  patented  Efratom  invention)  is 
used  in  the  Model  FRK.  Instead  of  using  a separate  temperature  sensor, 
the  current  that  is  used  to  heat  the  microwave  cavity  and  resonance  cell 
is  sensed  (this  current  decreases  as  ambient  temperature  increases)  and 
a signal  derived  from  this  is  used  to  correct  the  output  frequency  by 
changing  the  C-field  current.  In  fact,  this  very  technique  was  used  when 
making  the  long-term  stability  measurements  described  in  the  previous  sec- 
tion. For  those  measurements,  the  overall  temperature  coefficient  of  the 
unit,  with  electronic  temperature  compensation,  was  roughly  5 x 10" id/°C. 
There  are,  however,  two  limitations  to  this  technique.  First,  the  improve- 
ment in  temperature  coefficient  is  limited  to  about  a factor  of  ten  in  ac- 
tual practice.  Second,  the  transient  response  of  the  output  frequency  for 
rapid  temperature  changes  is  not  much  improved  because  the  different  ther- 
mal time  constants  associated  with  different  parts  of  the  apparatus  re- 
main the  same.  Clearly,  it  is  more  desirable  to  reduce  the  uncompensated 
temperature  sensitivity  than  to  try  to  compensate  for  it  electronically. 

In  summary,  it  appears  possible  to  reduce  the  sensitivity  of  the  out- 
put frequency  to  changes  in  ambient  temperature  by  using  the  methods  des- 
cribed above  (which  are  relatively  straightforward).  This  should  allow 
an  improvement  of  about  an  order  of  magnitude  in  the  overall  device  temp- 
erature coefficient. 
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FUNDAMENTAL  DEVICE  IMPROVEMENT 


As  noted  in  the  introduction,  there  are  two  parallel  approaches 
that  have  been  taken  in  this  work.  The  first  is  concerned  with  frequency 
stability  measurements  and  the  empirical  improvement  of  device  stability 
without  regard  to  the  physics  of  device  operation;  this  approach  is  des- 
cribed in  the  previous  section.  The  second  approach  is  an  attempt  to 
characterize  and  improve  device  operation  as  a frequency  standard  by  gain- 
ing a fundamental  understanding  of  device  behavior.  This  latter  approach 
is  considered  here  (see  Figure  3). 

\ 

Areas  of  Concern 

There  are  two  main  areas  of  concern  that  have  important  implications  \ 

for  operation  of  the  pulsed  pumping  apparatus  as  a frequency  standard.  I 

They  are  the  existence  of  possible  residual  light  shifts,  and  the  depen- 
dence of  the  output  frequency  on  the  pulse  parameters  (and  other  experi-  < 

mental  variables  as  well).  These  areas  of  concern  are  discussed  in  the 

next  two  sections.  . 


Sensitivity  to  Pulse  Parameters 


Figures  10  to  12  show  how  the  output  frequency  of  the  pulsed  pumping 
apparatus  depends  on  various  pulse  parameters.  In  all  three  figures,  the 
relevant  curve  (labeled  Af/f)  is  for  closed-loop  operation  as  a frequency 
standard.  Since  the  absolute  frequency  is  not  of  interest  here,  only  rela- 
tive frequency  differences  are  presented.  As  an  example.  Figure  10  shows 
that  if  the  light  pulse  duration  is  increased  from  1.0  msec  to  2.5  msec, 
the  output  frequency  increases  fractionally  by  3.0  x 10" 10.  Experimental 
conditions  for  the  three  Af/f  curves  are  given  in  Table  3.  The  Rb  lamp 
and  cell  temperatures  were  the  same  in  all  three  cases. 
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FM  SIGNAL,  SFM  (ARB.  UNITS) 


SIGNAL  CORRECTED  FOR  DUTY  CYCLE, AA™ (ARB. UNITS) 


FM  SIGNAL,  SFM  ARB.  UNITS 


Figure  11.  FM  signal  and  relative  fractional  frequency  offset  as  a func- 
tion of  dark  interval.  See  Table  3 for  experimental  conditions. 
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MICROWAVE  PULSE  DURATION , T,,  (msec) 


Figure  12.  Relative  fractional  frequency  offset  as  a function  of  micro- 
wave  pulse  duration.  See  Table  3 for  experimental  conditions. 
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TABLE  3.  (CONTINUED) 


The  Af/f  curves  in  Figures  10  to  12  all  exhibit  a common  behavior. 
There  is  a region  where  the  frequency  is  changing  very  rapidly  with  the 
pulse  parameter,  and  another  region  where  the  change  in  frequency  is  much 
more  gradual.  This  type  of  behavior  seems  to  hold  generally,  and  is  not 
dependent  on  a particular  parameter  having  a specific  value.  It  is  there- 
fore possible  to  choose  the  operating  conditions  so  that  frequency  changes 
due  to  changes  in  pulse  parameters  are  minimized,  and  this  is  usually  done. 

This  procedure  is  probably  adequate  to  ensure  long-term  stabilities  of  sev- 
1 ? 

eral  parts  in  10  or  better  over  days  and  weeks,  but  if  stabilities  of 
13  14 

parts  in  10  and  10  are  to  be  obtained  over  the  same  periods  of  time, 
and  longer,  then  it  will  probably  be  necessary  to  use  a more  fundamental 
approach  which  seeks  to  reduce  the  sensitivity  of  the  output  frequency  to 
changes  in  the  pulse  parameters,  and  other  instrumental  parameters  as  well 
(see,  for  example,  the  discussion  of  the  "position-shift  effect"  in  the 
next  section). 

A question  of  some  importance  is,  do  the  Af/f  curves  result  from 
fundamental  physical  effects,  or  are  they  some  spurious  artifact  of  the  ex- 
perimental setup?  It  is  difficult  to  answer  this  question  completely,  but 
it  is  possible  to  eliminate  a large  class  of  possible  artifacts.  We  have 
done  this  by  operating  the  apparatus  open  loop  (the  loop  is  broken  between 
the  integrator  and  the  VCXO  in  Figure  8)  and  observing  the  dc  component  of 
the  signal  from  the  output  of  the  phase  detector.  (In  effect,  we  are  using 
• the  signal  processing  section  of  Figure  5.)  The  VCXO  is  voltage  tuned  (as 

in  Figure  5)  so  that  the  6.8  GHz  is  way  off  resonance,  and  pulse  parameters 
are  then  varied.  Under  these  conditions,  any  spurious  effects  should  show 
up  as  a detectable  change  in  the  dc  output  level  of  the  phase  detector.  No 
detectable  change  is  observed.  The  same  result  is  also  obtained  when  the 
6.8  GHz  is  absent  altogether.  The  absence  of  any  change  in  the  dc  level  as 
. the  pulse  parameters  are  varied,  eliminates  a large  class  of  spurious  app- 

aratus effects  (e.  g.,  ground  loop  and/or  rf  pickup  induced  signals  in  the 
signal  processing  electronics  that  originate  in  the  pulsing  circuitry). 

When  this  type  of  experiment  is  carried  out  with  the  6.3  GHz  tuned  to  the 
peak  of  the  Rb  resonance,  dc  level  shifts  are  observed,  but  in  this  case  it 
is  not  possible  to  distinguish  fundamental  effects  from  spurious  ones.  To 
summarize,  we  have  carried  out  tests  to  detect  spurious  effects,  such  as 
ground  loops  and  rf  pickup,  and  have  detected  none.  Therefore,  in  so  far 
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as  we  can  tell,  the  variation  of  the  output  frequency  as  a function  of 
the  pulse  parameters  is  a real  effect,  and  is  not  just  an  electronic 
artifact  of  the  apparatus. 

Position-Shift  Effect 

The  so-called  "position-shift  effect"  has  been  studied^  by  Risley  and 
Busca.  This  effect  is  due  to  inhomogeneous  broadening  in  the  physics  pack- 
age, combined  with  nonuniform  excitation  (optical  pumping)  and  interroga- 
tion of  the  rubidium  atoms  as  discussed  below. 

There  are  two  mechanisms  that  can  cause  inhomogeneous  broadening  of 
the  rubidium  resonance  signal  from  the  physics  package.  The  first  is  in- 

homgeneity  of  the  static  magnetic  field.  Since  the  static  field  varies 

87 

over  the  volume  of  the  Rb  resonance  cell,  Rb  atoms  in  different  parts  of 
the  cell  will  have  slightly  different  resonant  frequencies  (Zeeman  effect). 
Another  effect  which  produces  a similar  end  result  is  the  light-shift  ef- 
fect. Due  to  the  light-filtering  action  of  the  Rb  atoms  in  the  resonance 
cell,  the  intensity  and  spectral  distribution  of  the  pumping  light  varies 

throughout  the  resonance  cell.  In  a cw  device  this  gives  rise  to  a light 

87 

shift  that  varies  throughout  the  resonance  cell.  That  is,  Rb  atoms  at 

different  locations  in  the  cell  will  have  different  resonant  frequencies 

due  to  the  light-shift  alone.  In  summary,  there  is  a nonuniform  distribu- 
87 

tion  of  Rb  atom  resonant  frequencies  throughout  the  resonance  cell  that 
is  caused  by  two  effects.  The  most  important  in  our  unit  is  due  to  the  in- 
homogeneity of  the  static  magnetic  field.  The  other,  which  should  be  ab- 
sent (or  else  very  small)  in  a pulsed  pumping  apparatus,  is  due  to  the 
• light-shift  effect. 

Inhomogeneity  of  the  above  type  is  not  sufficient,  by  Itself,  to  pro- 
duce frequency  shifts.  However,  in  combination  with  inhomogeneity  in  the 
87 

processes  of  Rb  atom  excitation  and  interrogation,  frequency  shifts  can 

and  will  occur. ^ This  can  easily  be  seen  as  follows.  As  discussed  above, 
87 

Rb  atoms  in  different  regions  of  the  resonance  cell  have  different  re- 
sonant frequencies.  Also,  because  the  optical  pumping  (excitation  to  the 
F = 2 hyperfine  level  of  the  ground  electronic  state)  and  microwave  inter- 
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rogation  is  nonuniform  over  the  cell,  the  relative  contributions  of 

various  portions  of  the  resonance  cell  is  not  uniform  over  the  cell. 

The  observed  atomic  line  shape  is  therefore  a weighted  superposition 
87 

of  individual  Rb  line  shapes  that  are  slightly  shifted  in  frequency 

relative  to  each  other.  The  center  frequency  and  shape  of  the  observ- 

ed  (composite)  Rb  line  therefore  depends  on  the  details  of  all  these 

inhomogeneities.  If  any  of  these  inhomogeneities  is  subject  to  change 

in  any  way,  then  there  will  be  a corresponding  change  in  the  shape  and 

87 

the  frequency  at  the  peak  of  the  observed  (composite)  Rb  resonance 
line.  This  will,  generally  speaking,  result  in  a frequency  shift  of  the 
device  when  operated  as  a frequency  standard.  Even  if  the  peak  fre- 
quency of  the  observed  resonance  line  is  not  shifted,  changes  in  its 
shape  that  produce  line  asymmetry  can  still  cause  output  frequency 

shifts.  This  can  happen  because  the  servo  system  locks  the  VCXO  to  the 
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center  of  gravity  of  the  observed  Rb  resonance  line  so  that  the  shape 
is  therefore  important  also. 

For  our  apparatus,  the  line  width  (FWHM)  is  typically  about  550  Hz, 
-8 

or  about  8 x 10  of  the  rubidium  resonance  frequency,  whereas  the  spread 
87 

of  Rb  resonant  frequencies  throughout  the  resonance  cell  is  of  order, 
g 

parts  in  10  . These  effects  therefore  constitute  a small  perturbation  on 
the  observed  line  shape,  but  they  are  quite  large  compared  to  the  preci- 
sion with  which  the  frequency  of  the  device  is  to  be  stabilized  (~  1 x 
10  ^ and  better). 

An  example  of  the  type  of  frequency  shift  that  can  occur  due  to  the 
position-shift  effect  is  that  due  to  microwave  power.  When  the  microwave 
power  applied  to  the  microwave  cavity  changes,  then  some  regions  of  the 
resonance  cell  will  be  interrogated  more  strongly,  and  other  regions  more 
weakly.  If  the  regions  that  are  interrogated  more  strongly  have  differ- 
ent resonance  frequencies  than  those  that  are  interrogated  more  weakly,  a 
frequency  shift  will  result. 

Risley  and  Busca^  have  presented  evidence  which  shows  that  the  posi- 
tion-shift effect  may  be  a limiting  factor  in  the  long-term  stability  of 
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j cw,  passive  Rb  frequency  standards.  In  a pulsed  optical  pumping  device, 

it  should  be  possible  to  eliminate,  or  reduce  to  a negligible  value,  fre- 
f quency  shifts  that  are  due  to  the  position-shift  effect.  This  can  be  done 

by  making  the  static  magnetic  field  more  uniform.  In  the  FRK  physics 
package  that  we  are  using  in  the  pulsed  pumping  apparatus,  this  field 
is  quite  nonuniform  and  a considerable  improvement  should  be  easy  to 
obtain.  In  principle,  if  the  static  magnetic  field  in  a pulsed  pumping 
apparatus  is  completely  homogeneous,  the  position-shift  effect  should  be 
entirely  eliminated.  This  is  not  true  for  a cw  apparatus  because  there 
is  still  inhomogeneous  broadening  due  to  light  shifts.  This,  therefore  is 
a definite  advantage  of  the  pulsed  pumping  method. 

Residual  Light  Shift 

A question  of  great  interest  from  a practical  point  of  view  (and 
otherwise)  is,  is  there  any  light  shift  when  pulsed  optical  pumping  is 
used?  For  the  apparatus  that  we  have  constructed,  there  is  a small  resi- 
dual light  shift.  On  the  other  hand,  there  is  evidence  which  strongly 
indicates  that  this  is  not  a true  light-shift  effect,  but  arises  instead 
from  a different  physical  effect.  For  this  reason  we  will  refer  to  this 
phenomenon  as  the  "pseudo  light-shift  effect"  (see  Figure  3).  This  resi- 
dual light  shift  is  discussed  below. 

An  experiment  was  conducted  using  two  different  rubidium  lamps,  #210 
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and  #110.  Lamp  #210  contains  pure  Rb,  while  lamp  #110  contains  both 
rubidium  isotopes.  The  spectral  distributions  of  the  light  from  these  two 
lamps  are  expected  to  differ  greatly  and  therefore  to  produce  significantly 
different  light  shifts  in  a gas  cell  rubidium  device  operating  in  the  cw 
mode.  For  operation  in  the  pulsed  mode  (with  no  light-microwave  overlap), 
the  light  shifts  are  expected  to  be  either  absent,  or  else  to  be  very  much 
smaller  in  magnitude  than  the  light  shifts  for  cw  operation  (but  of  the 
same  sign  in  each  case). 

The  experiment  consisted  of  first  placing  lamp  #210  in  the  apparatus. 
For  operation  in  the  cw  mode,  the  fractional  frequency  of  the  apparatus. 
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relative  to  a fixed  reference,  was  -2.6  x 10  . Lamp  #210  was  then  re- 

moved and  replaced  by  lamp  #110.  With  the  apparatus  still  operating  cw, 

the  fractional  frequency,  relative  to  the  same  fixed  reference,  was  +3.9  x 
_q 

10  for  lamp  #110.  The  total  frequency  shift  in  the  cw  mode,  upon  chang- 

_g 

ing  lamps,  was  therefore  +6.5  x 10  . When  the  apparatus  was  operated  in 

the  pulsed  mode,  under  the  same  conditions  with  the  same  fixed  reference, 
the  fractional  frequency  offsets-  were  -2  x 10"^  for  lamp  #210,  and  -4  x 
10'10  for  lamp  #110.  The  total  frequency  shift  in  the  pulsed  mode,  upon 
changing  lamps,  was  therefore  -2  x 10"^.  Similar  experiments  with  the 
lamp  pulsed  and  the  microwave  operated  cw  (with  approximately  50  % light- 
microwave  overlap)  gave  frequency  offsets  and  shifts  that  lay  - halfway 
between  the  values  for  the  fully  pulsed  and  the  fully  cw  modes. 

The  results  of  this  experiment  are  consistent  with  the  hypothesis 

_9 

that  the  frequency  shift  of  +6.5  x 10  in  the  cw  mode  is  due  mostly  to  a 

change  in  the  light  shift  when  lamp  #210  is  replaced  by  lamp  #110.  If 

this  is  a valid  hypothesis,  then  the  shift  of  -2  x 10-1^  in  the  pulsed 

mode  (under  the  same  conditions),  must  be  interpreted  as  due  to  some  ef- 

. 15 

feet  other  than  a true  light  shift  because  of  the  opposite  sign. 

A different  way  to  test  for  light  shifts  in  the  cw  and  pulsed  modes 
1 is  to  measure  the  frequency  shift  in  each  mode  as  a function  of  light  in- 

tensity for  a given  lamp  (by  inserting  light  attenuating  filters  in  the 
light  beam).  Such  measurements  have  been  made  in  the  pulsed  mode  for  two 
different  light  pulse  durations,  but  no  such  measurements  have  been  made 
in  the  lw  mode  for  the  same  experimental  conditions  (e.g.,  lamp  tempera- 
ture and  cell  temperature).  The  results  for  the  pulsed  mode  are  shown  in 
Figure  13,  where  it  is  seen  that  the  fractional  frequency  shift  is  a lin- 
ear function  of  light  intensity.  Although  two  straight  lines  are  shown  in 
the  figure,  the  uncertainties  in  the  plotted  points  and  circles  are  of  the 
order  of  the  difference  between  the  two  lines.  We  must  therefore  conclude 
that  the  shift  is  sensibly  independent  of  light-pulse  duration. 

Additional  measurements  were  made  with  each  lamp  in  the  pulsed  pump- 
ing mode  by  inserting  a light  attenuating  filter  in  the  light  beam  at  the 
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RELATIVE  LIGHT 
INTENSITY  CHANGE,  I/I0 


Figure  13.  Pseudo  light-shift  dependence  on  light  intensity.  Open  cir- 
cles for  T.  = 2.0  ms,  dots  for  Ti  = 1.1  ms;  other  parameters  are  given 

in  Table  3. 
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point  where  it  leaves  the  Rb  lamp.  It  has  been  previously  established 
that  this  filter  attenuates  the  Rb  and  D2  lines  by  22  %.  The  effect 
of  the  filter  is  to  produce  a shift  in  the  output  frequency  of  the  unit 
of  -1.6  x 10"10.  The  same  result  is  obtained  for  both  lamps  (#210  and 
#110),  whereas  for  a true  light  shift  we  expect  the  shift  for  lamp  #210 
to  be  about  half  that  for  lamp  #110,  and  of  the  opposite  sign. 

As  another  check  on  whether  this  is  a true  light  shift,  we  can  exam- 
ine how  this  filter  shift  varies  with  the  duration  of  the  dark  period 

(T,).  The  theory  (with  supporting  experimental  data)  developed  by  Ba- 
u 15 

zarov  and  coworkers  predicts  that  a true  residual  light  shift  decreases 

in  magnitude  as  T^  is  increased.  Therefore,  we  measured  the  filter  shift 

for  Td  = 1.8  msec  and  T^  = 2.8  msec  (other  parameters  were  T^  = 1.1  msec, 

T = 0.14  msec,  and  T = 1.5  msec).  It  was  found  that  the  filter  shift 
S M -in 

was  the  same,  and  equal  to  -1.6  x 10  in  both  cases.  This  is  contrary 

to  the  behavior  expected  for  a true  light  shift,  as  noted  above. 

On  the  basis  of  these  experiments,  it  seems  highly  unlikely  that 
this  is  a true  light  shift.  At  the  present  time,  the  origin  of  this 
pseudo  light  shift  is  not  understood.  One  possibility  is  that  it  is  due, 
at  least  in  part,  to  the  position-shift  effect:  as  the  light  intensity  is 
changed  the  spatial  distribution  of  excitation  (optical  pumping)  in  the 
resonance  cell  is  also  changed,  and  there  is  a resulting  frequency  shift 
(that  is  really  due  to  the  inhomogeneous  broadening  by  the  static  magne- 
tic field  and  has  nothing  whatsoever  to  do  with  a true  light  shift).  One 
way  of  checking  this  hypothesis  is  to  repeat  the  above  measurements  with  a 
homogeneous  static  magnetic  field.  If  the  pseudo  light-shift  effect  per- 
sists unchanged,  then  the  position-shift  effect  is  not  the  cause. 

In  summary,  there  is  a small  residual  "light  shift"  when  our  pulsed 
optical  pumping  apparatus  is  operated  as  a frequency  standard.  Various 
experiments  have  shown  that  this  shift  does  not  behave  like  a true  light 
shift.  For  this  reason,  we  have  called  it  the  "pseudo  light  shift."  Al- 
though the  origin  of  the  pseudo  light  shift  is  not  presently  known,  it  is 
possible  that  it  may  arise  from  the  position-shift  effect.  A way  of  test- 
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ing  this  hypothesis  has  been  presented.  The  dependence  of  the  pseudo 

light  shift  on  light  intensity  has  been  measured  (Figure  13),  and  from 

these  data  it  is  estimated  that  a 1 % change  in  light  intensity  will  pro- 

-12 

duce  a frequency  shift  of  6 x 10  . From  this,  it  is  clear  that  if  the 

long-term  stability  of  the  pulsed  pumping  apparatus  is  to  be  improved  be- 

12 

yond  its  present  value  of  parts  in  10  , the  pseudo  light  shift  will  have 

to  be  reduced  by  at  least  an  order  of  magnitude.  This  should  be  possible 
if  the  pseudo  light  shift  is  a manifestation  of  the  position-shift  effect, 
but  further  work  will  be  required  to  establish  this. 

Understanding  Device  Operation 

The  previous  section  (Areas  of  Concern)  has  shown  that  there  is  a 
definite  need  to  obtain  a fundamental  understanding  of  device  operation. 
Reference  to  Figure  3 shows  that  some  progress  along  this  line  has  already 
been  made.  This  is  discussed  below. 

Theoretical  Model  of  Device  Operation 

The  development  of  a theoretical  model  of  device  operation  must  pro- 
ceed in  two  steps.  The  first  is  the  prediction  of  how  the  observed  rubi- 
dium resonance  signal  depends  on  the  parameters  of  the  apparatus,  and  the 
second  is  the  prediction  of  how  the  output  frequency  depends  on  the  same 
parameters. 

15 

The  theory  of  Bazarov  & coworkers  treats  both  aspects  of  device 
operation  with  emphasis  on  the  latter,  but  it  also  neglects  an  important 
aspect  of  device  operation  that  is  important  for  our  apparatus,  namely, 
inhomogeneous  broadening  (such  as  occurs  in  the  position-shift  effect). 

For  this  reason,  we  have  chosen  to  develop  our  own  theoretical  treatment, 
as  described  below.  So  far,  we  have  not  attempted  to  formulate  a theory 
of  the  line  shape,  but  our  theory  could  be  extended  to  include  this. 
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Starting  with  the  basic  quantum  mechanical  equations  that  describe 
the  optical  pumping  process,  we  have  developed  a theoretical  model  that 
predicts  how  the  observed  signal  depends  on  various  apparatus  parameters 
other  than  microwave  frequency,  especially  the  pulse  parameters.  This 
theory  is  briefly  outlined  in  Appendix  A.  In  the  next  section  (AM  meas- 
urements), the  results  of  an  experimental  investigation  of  device  proper- 
ties are  presented,  and  compared  with  the  predictions  of  this  theory.  Good 
agreement  between  theory  and  experiment  is  obtaned. 

AM  Measurements 


The  AM  measurements  fall  into  two  categories,  as  can  be  seen  in  Fig- 
ure 3.  Each  category  is  discussed  separately  below. 

Line  Shapes 

87 

Observed  Rb  resonance  (AM)  line  shapes  are  shown  in  Figure  14  for  a 

variety  of  experimental  conditions  (as  itemized  in  Table  3).  The  lamp  and 

cell  temperatures  are  identical  for  all  curves  except  curve  D.  For  this 

reason  we  will  discuss  curve  D last.  All  curves,  except  curves  A and  F, 

were  obtained  using  lamps  that  contained  both  isotopes  of  Rb.  Curves  A and 

87 

F were  taken  using  a lamp  containing  only  Rb.  The  type  of  lamp  used  is 
not  expected  to  have  much  effect  on  the  line  shape  on  the  scale  shown  in 
Figure  14. 

Curve  B shows  the  line  shape  for  cw  operation.  In  this  case  the  line 
28 

width  is  approximately  1100  Hz.  However,  when  the  apparatus  is  operated 
in  the  pulsed  mode  it  is  possible  to  obtain  line  widths  as  small  as  360  Hz 
(curve  F).  This  shows  that  under  cw  conditions  light  relaxation  makes  a 
significant  contribution  to  the  resonance  line  width.  This  also  demon- 
strates that  in  the  pulsed  mode  the  microwave  radiation  really  does  inter- 
rogate the  atoms  when  they  are  in  the  dark.  From  the  uncertainty  principle, 
it  is  expected  that  the  line  width  in  the  pulsed  mode  should  vary  as  1/T^, 
where  T^  is  the  duration  of  the  exciting  microwave  radiation.  For  line 
widths  larger  than  360  Hz  this  is  the  case,  and  the  data  show  that  the  line 
width  ~ 0.9/T^.  In  this  regard,  it  is  interesting  to  note  that  when  T^  be- 
comes very  small,  the  line  width  becomes  quite  large.  This  is  demonstrated 
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Figure  14.  Observed  Rb  resonance  AM  line  shapes.  All  curves  are  normal- 
ized to  the  same  value  at  Af/f  = 0.  Experimental  conditions  are  given 

in  Table  3. 
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by  curve  A where  the  line  width  for  pulsed  operation  with  = 0.55  msec 
is  about  1700  Hz,  which  is  greater  than  the  cw  line  width  (curve  B). 

The  lower  limit  on  the  line  width  in  the  pulsed  mode  is  set  by  in- 
homogeneity broadening  of  the  static  magnetic  field,  and  various  relax- 
ation processes.  It  is  estimated  that  the  magnetic  field  inhomogeneity 
contributes  less  than  50  Hz.  The  contribution  due  to  relaxation  is  larger 
than  this,  and  can  be  estimated  from  measurements  (see  the  next  section) 
of  the  relaxation  time  for  the  decay  of  the  F=2,  F=1  population  difference 
(Tj  in  NMR  jargon).  For  a cell  temperature  of  75  °C,  T^  = 1.2  msec. 

Since  we  do  not,  at  present,  have  a theory  of  the  line  shape,  it  is  not 
clear  whether  T^  or  T?  (relaxation  time  for  decay  of  the  oscillating  mag- 
netization) or  some  combination  of  the  two,  is  appropriate  to  use  to  com- 
pute the  relaxation  contribution  to  the  line  width.  However,  the  differ- 
ences are  expected  to  be  small  because  usually  T?  «T^.  If  we  use  the  ex- 
perimental value  of  T^,  the  contribution  due  to  relaxation  is  260  Hz.  To 
this  must  be  added  (directly)  the  contribution  due  to  inhomogeniety  broad- 
ening. This  gives  a final  estimate  for  the  total  line  width  of  approxi- 
mately 300  Hz  which  is  in  good  agreement  with  the  experimental  observa- 
tions (experimental  lower  limit  <360  Hz). 

We  now  discuss  curve  0.  For  it,  two  sets  of  data  were  taken,  differ- 
ing only  in  the  values  of  T^.  The  dots  are  for  T?  = 1.8  msec,  and  the 
open  circles  for  T^  = 1.1  msec.  It  is  not  known  if  the  slight  difference 
observed  between  these  two  sets  of  points  is  significant.  The  wing  of  the 
line  for  curve  D is  also  shown  separately  in  Figure  14.  The  subsidiary 
maximum  is  a real  effect  and  has  been  observed  in  other  line  shapes  (not 
shown  in  Figure  14)  when  the  apparatus  is  operated  in  the  pulsed  mode. 

This  effect  is  not  present  when  the  apparatus  is  operated  in  the  cw  mode. 

Subsidiary  maxima  also  occur  in  molecular  beam  experiments  when  the  oeam 
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is  monochromatic  in  velocity.  ’ In  both  cases  the  basic  effect  is  due 

30 

to  the  microwave  (or  radiofrequency)  radiation  being  applied  as  a pulse. 

It  might  be  thought  that  it  would  be  advantageous  to  operate  the 
pulsed  pumping  apparatus  using  the  narrowest  possible ■ line  width.  This  is 
not  necessarily  the  case,  however,  because  the  signal  strength  decreases 
approximately  exponentially  with  Ty,  whereas  the  line  width  decreases  only 
as  mv 
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AM  Signal  as  a Function  of  Pulse  Parameters 

The  AM  signal  has  been  studied  as  a function  of  the  pulse  parameters 
T£,  Td  and  The  variation  as  a function  of  light  pulse  duration  is 

shown  in  Figure  15.  The  dots  are  the  experimental  data  points,  and  the 

AM 

curve  labeled  SAM  passes  through  these  points.  The  curve  labeled  AA 

is  obtained  (except  for  the  overall  normalization  which  is  arbitrary)  by 

multiplying  the  function  SflM(T0)  by  (T0  + T.).  According  to  Eq.  (A24), 

mm  x,  x,  a ^ 

this  gives  the  AM  signal  corrected  for  duty  cycle,  AA  . 

AM 

The  behavior  of  AA  in  Figure  15  is  as  expected  on  the  basis  of  the 

theory  developed  in  Appendix  A.  Starting  with  a narrow  light  pulse  the 

AM  signal  AA^  (corrected  for  duty  cycle  effect)  increases  rapidly  at 

first  with  light  pulse  duration,  and  subsequently  approaches  saturation. 

Physically,  this  can  be  understood  as  follows.  If  the  lamp  is  suddenly 

87 

turned  on  (the  microwave  is  assumed  to  be  off  for  the  moment),  Rb  atoms 
will  be  pumped  out  of  the  F=1  hyperfine  level  and  into  the  F=2  level.  This 
causes  the  population  of  the  F=1  level  to  decrease  exponentially  with  time, 
approaching  an  asymptotic  value  that  depends  on  the  intensity  of  the  light 
and  the  relaxation  rate.  Since  the  light  absorption  is  proportional  to  the 
population  of  the  F=1  level,  the  light  absorption  will  be  greatest  a*  first 
and  will  also  decrease  exponentially  with  pumping  time.  This  type  of  be- 
i havior  is  shown  schematically  in  Figure  16  for  two  cases.  In  the  first 

case,  the  light  pulse  duration  is  T^  = T.  During  one  half  of  the  10  Hz 
modulation  cycle,  the  microwave  radiation  is  not  present  between  light 
pulses.  In  the  steady  state,  this  gives  the  lower  curve  shown  in  Figure 
16.  During  the  other  half  of  the  10  Hz  modulation,  the  microwave  pulses 

5 

are  present  between  the  light  pulses.  Under  steady  state  conditions,  this 

gives  the  upper  curve  of  Figure  16.  In  this  case,  the  AM  signal  corrected 
AM 

for  duty  cycle,  AA  , is  proportional  to  the  shaded  area  that  is  between 
the  two  curves  in  the  interval  T^  = 0 to  T^  = T.  If  the  light  pulse  dura- 
tion is  increased  from  T = T to  T = T'then  this  adds  the  additional  shad- 

ed  area-31  between  T and  l'.  Looking  at  Figure  16,  it  is  now  clear  why  the 
AM 

AM  signal  (AA  in  Figure  15)  increases  rapidly  at  first  and  then  approach- 
es saturation  as  T^  is  increased. 

The  behavior  cf  the  AM  signal  as  a function  of  the  dark  time  Td  has 
also  been  investigated  experimentally.  All  experimental  parameters  except 
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Figure  15.  Measured  AM  signal  as  a function  of  light  pulse  duration. 
See  Table  3 for  experimental  conditions. 
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AM  SIGN 


T,  are  held  fixed.  Also,  T = 2.0  msec  which  ensures  that  we  are  oper- 

o *■  am 

ating  in  the  saturation  region  for  AA  (see  Figure  15),  and  Td  is  varied 

in  such  a way  that  there  is  no  light  - microwave  overlap  (T^  > T$  + T ) . 

Under  these  conditions,  Eq.  (A23)  applies.  This  equation  shows  that  in 

AM 

the  light-pulse-duration  saturation  limit,  the  AM  signal  AA  (corrected 

-Y  T 

for  duty  cycle)  is  proportional  to  e 1 d.  This  allows  a determination 

of  the  relaxation  time  T^  = l/y^  which  is  the  time  constant  for  the  decay 

of  the  population  difference.  Figure  17  shows  the  experimental  results 

for  a change  in  AM  signal  of  > 10:1.  According  to  Eq.  (A23),  the  curve 

of  Figure  17  should  be  a single  straight  line.  Instead,  the  data  are  fit 

by  two  straight  lines,  A-A  and  B-B.  For  A-A,  T^  = 1.67  msec,  whereas  for 

B-B,  T,  = 1.82  msec.  While  this  effect  may  be  an  experimental  artifact, 

1 32 

it  is  also  possible  that  it  is  real.  This  is  further  discussed  below. 

Table  4 summarizes  both  measured  and  calculated  values  of  Tj  for  var- 
ious resonance  cell  temperatures.  The  resonance  cell  temperature  stated 
in  the  table  as  66°C  is  estimated  to  be  accurate  to  better  than  ±0.5  °C 
at  the  location  of  the  thermocouple.  However,  there  is  a qradient,  per- 
haps as  large  as  1 °C,  across  the  resonance  cell.  The  remaining  two  val- 
ues of  temperature  in  Table  4 were  measured  very  early  in  this  work  and 

have  large  uncertainties  (~  +2  °C). 

87  87 

The  contribution  of  Rb-  Rb  spin-exchange  relaxation  to  the  total 
Tj  is  given  in  column  three  of  Table  4.  These  values  were  obtained  by 
calculating  the  total  Rb  pressure  in  the  resonance  cell  from  the  measured 
resonance  cell  temperature  using  the  equation,  log^g  P(torr)  = -4560/Tc 
(°K)  + 12.00  - 1.45  log-0  Tc(°K).  The  spin-exchange  contribution  was  then 
calculated  from  the  experimental  result  of  reference  33,  namely,  1/CTj)^ 

= 6.75  x 107  sec"1  x P(torr). 

In  principle,  the  total  T^  cannot  exceed  that  due  to  spin- exchange 
alone.  For  our  apparatus  and  operating  conditions,  it  expected  tnat  con- 
tributions to  the  total  y^  due  to  buffer  gas  relaxation  and  wall  relaxa- 
tion will  be  quite  small.  That  is,  our  experimental  values  of  T^  should 
be  very  nearly  equal  to  the  value  of  Tj  due  to  spin-exchange  alone.  The 
discrepancies  that  occur  in  Table  4 may  be  attributed  to  numerous  system- 
atic errors  that  can  plague  this  type  of  measurement:  especially  errors  in 
measuring  temperature,  and  the  existence  of  temperature  gradients  (already 
discussed  above);  the  fact,  previously  reported,  that  Tj  depends  on  the 
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AM  SIGNAL  CORRECTED  FOR  DUTY  CYCLE, A Amm(arb. units) 


0 12  3 4 5 6 

DARK  TIME,  Tj  (msec.} 


Figure  17.  Measured  AM  signal  as  a function  of  dark  time.  See  Table  3 for 

experimental  conditions. 
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TABLE  4.  MEASURED  AND  CALCULATED 
RELAXATION  TIMES  (T^. 


Resonance  cell  Relaxation  time,  T.  (msec) 

07  * 07 

temperature,  Tc  total  Rb-  Rb  spin-exchange 

(°C)  experimental3  calculated 


66 

1.67c,d 

1.95 

1.82c,e 

68 

1 . 74f 

1.64 

75 

1.199 

0.91 

aPresent  work 

bCalculated  from  the  experimental  data  of  reference  33. 

cData  taken  for  11:1  change  in  AAAM,  and  could  not  be 
fit  by  a single  straight  line. 

dFrom  slope  of  straight  line  A-A,  Figure  17. 

eFrom  slope  of  straight  line  B-B,  Figure  17. 

f AM 

Data  taken  for  3:1  change  in  AA  , and  could  be  fit 

by  a single  straight  line. 

gData  taken  for  5:1  change  in  AAAM,  and  could  be  fit 
by  a single  straight  line. 


history  of  the  resonance  cell  and  the  location  of  the  rubidium  within 
it,  in  addition  to  the  cell  temperature;  the  same  considerations  just 
mentioned,  applied  to  the  experimental  data  and  analysis  of  reference  33. 

One  other  possible  source  of  discrepancy  in  comparing  columns  two 

and  three  of  Table  4 should  also  be  mentioned.  The  measurements  of  ref- 
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erence  33  were  made  using  resonance  cells  filled  with  Rb.  Our  reson- 
ance cell  contains  both  isotopes  of  Rb,  so  that  spin-exchange  relaxation 
of  the  87Rb  is  due  to  87Rb-85Rb  collisions  as  well  as  to  87Rb-87Rb  col- 
lisions. Since  the  nuclear  spin  can  have  a significant  effect  on  relax- 
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ation  due  to  spin  exchange,  it  is  possible  that  there  are  two  relaxa- 

87  87  87  85 

tion  times,  one  for  Rb-  Rb  collisions  and  another  for  Rb-  Rb  col- 
lisions. If  this  is  the  case,  then  it  would  not  be  surprizing  to  observe 
two  different  decay  curves.  This  might  be  the  origin  of  the  two  decay 
curves  seen  in  Figure  17. 

The  AM  signal  SAM  as  a function  of  microwave  pulse  duration  Ty  is 
shown  in  Figure  18  as  a solid  curve.  (In  this  case,  no  correction  for 
duty  cycle  is  required  because  the  duty  cycle  factor  (T^  + T^)-1  is  con- 
stant.) An  intriguing  characteristic  of  this  curve  is  the  presence  of  two 
maxima,  and  because  of  this  it  is  of  interest  to  compare  the  experimental 
result  with  theoretical  predictions.  Since  the  experimental  conditions 
(Table  3)  were  such  that  the  light  pulse  duration  was  in  the  saturation 
region,  the  simplified  theory  in  Appendix  A can  be  used.  Three  levels  of 
sophistication  are  possible.  We  discuss  them  in  order  of  increasing  com- 
plexity. 

First,  we  consider  the  result  obtained  in  the  very  simple  case  where 
the  microwave  magnetic  field,  and  the  optical  pumping  are  uniform  through- 
out the  resonance  cell.  The  predicted  dependence  of  AM  signal  S^M  is  giv- 
en by  Eq.  (A21)  multiplied  by  the  factor  1/(T^  + Td),  as  per  Eq.  (A24). 
This  equation  is  not  plotted  in  Figure  18  because  it  is  obvious  by  inspec- 
tion that  it  gives  a poor  fit  to  the  data.  It  does,  however  predict  suc- 
cessive maxima,  but  the  minima  are  predicted  to  be  much  too  deep,  and  suc- 
cessive peaks  are  too  far  apart. 

The  next  degree  of  sophistication  is  Eq.  (A27),  which  is  obtained  by 
averaging  Eq.  (A21)  over  the  resonance  cell,  as  explained  in  Appendix  A. 
Physically,  Eq.  (A27)  holds  for  a microwave  magnetic  field  that  varies  as 
half  of  a sine  wave  over  the  length  of  the  resonance  cell,  and  for  con- 
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MICROWAVE  PULSE  DURATION,  TM(msec) 


Figure  18.  AM  signal  as  a function  of  microwave  pulse  duration.  Solid 
curve  (and  points)  are  measured;  other  curves  are  theoretical.  See  Ta- 
ble 3 for  experimental  conditions. 
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stant  optical  pumping  over  the  length  of  the  cell.  This  equation, 
when  normalized  and  fit  to  the  first  peak  in  Figure  18  gives  the  dash- 

dot  curve  shown  there.  In  plotting  this  curve.  Ay  = 210  sec  * was  as-  - 

sumed.  This  value  of  Ay  is  the  value  that  results  if  it  is  assumed  that 

87  87  34 

all  relaxation  is  due  to  Rb-  Rb  spin-exchange  so  that  T£  = (8/5)  T^,  i 

and  the  experimental  value  of  T^  from  Table  4 is  used  (T^  = 1.8  msec). 

The  dash-dot  curve  is  an  improvement  over  the  prediction  of  Eq.  (A21),  but 
it  still  suffers  from  the  same  basic  drawbacks. 

The  final  degree  of  sophistication  that  we  have  considered' results 
from  assuming  that  the  optical  pumping  is  much  greater  at  the  end  of  the 
resonance  cell,  than  at  the  entrance7  (see  Figure  A3),  and  that  the  micro- 

wave  field  variation  is  the  same  as  before.  As  noted  in  Appendix  A,  numer-  ^ 

ical  calculation  is  required  under  these  assumptions.  The  result  for  Ay  = 

210  sec"*  is  shown  in  Figure  18  as  the  large-dash  curve.  For  Ay  = 0,  the 

small-dash  curve  is  obtained  instead.  The  best  fit  is  obtained  for  Ay  = ; 

0.  The  parameter  8.  is  adjusted  to  give  the  best  fit.  For  both  curves, 

- -1  0 

8q  = 1600  sec  was  used.  This  corresponds  to  a z-component  of  microwave  * 

magnetic  field  amplitude  at  the  center  of  the  cell  of  approximately  360  pG.  ( 

Figure  18  shows  that  the  theory  is  over  simplified,  but  it  also  shows 
that  it  contains  a great  deal  of  the  truth.  The  deficiencies  of  the  theory 
probably  lie  more  in  the  assumptions  made  when  averaging  over  the  resonance 
cell,  than  in  the  basic  (quantum  mechanical)  theory.  For  example,  since  \ 

the  static  magnetic  field  is  not  uniform  in  magnitude  or  direction  (as  as- 
sumed in  all  cases  above),  the  z-component  of  the  microwave  field  Hz  is  not 
the  only  component  that  is  effective  in  inducing  the  0-0  hyperfine  transi-  j 

tion.  Also,  Hz  depends  strongly  on  azimuthal  angle  as  well  as  on  z,  but  I 

the  former  was  ignored  in  the  theoretical  treatment.  It  is  likely  that  if 
these  effects  were  taken  into  account,  the  sharp  dips  in  the  theoretical 
curves  of  Figure  18  would  be  smoothed  out. 

I 

FM  Measurements 

I 

The  variation  of  the  FM  signal  SpM  was  also  investigated  as  a func-  , 

tion  of  the  pulse  parameters.  The  procedure  used  was  to  open  the  servo 
loop  between  the  integrator  and  the  VCX0  (see  Figure  8);  to  replace  the 
VCX0  by  a more  stable  10  MHz  source;  to  tune  this  source  slightly  off  re- 
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sonance,  thereby  giving  a nonzero  FM  signal  in  accordance  with  Figure  9 
(care  was  taken  to  remain  on  the  linear  portion  of  the  discriminator  char- 
acteristic of  Figure  9);  to  replace  the  integrator  of  Figure  8 by  the  low 
pass  filter  and  chart  recorder  of  Figure  5.  The  variation  of  Sp^  was  then 
studied  as  a function  of  light  pulse  duration  (Tf ) and  dark  time  (Td). 

The  results,  which  are  shown  in  Figures  10  and  11,  are  qualitatively  simi- 
lar to  the  AM  results.  In  Figure  10,  Sp^  increases  rapidly  at  first, 
reaches  a maximum,  and  then  falls  off  gradually.  AA™,  the  FM  signal  cor- 
rected for  duty  cycle,  increases  rapidly  at  first  and  then  reaches  a sat- 
uration condition.  This  is  very  similar  to  the  corresponding  AM  behavior 
of  Figure  15.  Quantitatively,  differences  are  observed.  For  example,  the 
peak  of  the  SpM  curve  is  at  Tp  = 2.0  msec,  whereas  the  peak  of  the  S^M 


curve  is  at  T 


1.7  msec. 


Figure  11  shows  the  FM  signal  SpM  as  a function  of  dar  ■>me  T^.  If 
this  curve  is  plotted  on  semilog  paper,  the  experimental  points  cannot  be 
precisely  fit  by  one  or  two  straight  lines,  but  an  approximate  fit  can  be 
obtained  using  a single  straight  line.  For  this  straight  line,  Tp  = 1.76 
msec.  This  is  very  nearly  equal  to  the  average  of  the  two  values  of  Tp 
obtained  from  the  AM  measurements  at  the  same  resonance  cell  temperature 
of  66°C  (see  Table  4). 

When  FM  is  used,  the  derivative  of  the  resonance  curve  is  obtained 
for  small  modulation  index.  However,  signal/noise  considerations  require 
that  a large  modulation  index  be  used,  in  which  case  the  FM  line  shape  is 

or  or 

not  simply  related  to  the  Rb  resonance  curve.  * For  this  reason,  it  is 
difficult  to  interpret  quantitatively  the  FM  measurements. 


SUMMARY  AND  CONCLUSIONS 


1.  Technical  Feasibility 

The  design  philosophy  adopted  in  the  present  work  is  that  it  is  not 
profitable  to  conduct  laboratory  experiments  using  principles  and  methods 
that  cannot  subsequently  be  applied  to  the  development  of  a compact,  light 
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weight,  low  power  frequency  standard.  In  keeping  with  this  design  philos- 
ophy, al 1 experimental  work  described  in  this  report  has  been  carried  out 
using  a very  small  physics  package  ( consisting  of  lamp,  "integrated  cell," 
microwave  cavity,  magnetic  shield,  photodetectors,  and  lamp  electronics) 
that  occupies  a total  volume  of  260  cc  (about  equal  to  the  volume  of  an  or- 
ange.) Using  this  physics  package,  it  has  been  demonstrated  that  the 
method  of  pulsed  optical  pumping  is  technically  feasible  in  the  sense  that 
signal -to-noise  ratios  can  be  obtained  that  are  adequate  for  operating  the 
device  as  a frequency  standard.  This  is  possible  because  of  a dual  photo- 
cell compensation  method  that  was  developed  in  the  course  of  this  inves- 
tigation. 

2.  Short-Term  Frequency  Stability  (x  = 1 to  100  sec. ) 

The  short-term  frequency  stability  has  been  measured  for  the  pulsed 
optical  pumping  frequency  standard  described  above.  The  result  is 

ay(x)  s 2 x 10’^  i”5,  1 < t < 100  sec. 

When  the  pulsed  pumping  apparatus  is  operated  cw  (no  pulsing)  the  improve- 
ment in  Oy  is  less  than  a factor  of  two. 

It  should  be  pointed  out  that  the  Oy  obtained  for  pulsed  pumping  is 
comparable  to  that  of  present-day  commercial  rubidiums,  and  exceeds  that 
of  commercial  cesiums  (see  Table  1)  for  the  same  averaging  times.  Also, 
it  is  possible  that  this  Oy  could  be  improved,  perhaps  by  a factor  of  two, 
with  further  work. 

3.  Long-Term  Frequency  Stability  (t  = 24  hours) 

Preliminary  measurements  of  long-term  stability  indicate  that  the  24 

hour  average  frequency  of  the  pulsed  optical  pumping  frequency  standard 

12 

does  not  change  by  more  than  several  parts  in  10  from  one  24  hour  period 
to  the  next,  and  may  actually  be  better  than  this.  At  the  present  time, 
our  verifiable  (by  measurements  traceable  to  NBS)  in-house  measurement  ca- 
pability is  < 4 x 10”12  for  t = 24  hours.  Please  see  NOTE  ADDED  IN  PROOF, 

p.  68. 


61 


4.  Temperature  Coefficient  of  Frequency 

Over  a 10  °C  range  of  ambient  temperature,  the  temperature  coeffi- 
cient (TC)  of  frequency  of  the  pulsed  optical  pumping  frequency  standard 
is  approximately  +5  x 10"^/°C.  This  is  within  a factor  of  two  of  the  TC 

3 

of  the  company's  commercial  rubidium.  For  the  pulsed  device,  both  the 
electronics  and  the  Rb  resonance  cell  are  temperature  sensitive,  and  it 
should  be  possible  to  make  significant  improvements  in  both  areas.  For 
example,  by  using  a different  combination  of  buffer  gases  in  the  resonance 
cell,  it  should  be  possible  to  reduce  the  resonance  cell  TC  by  nearly  a 
factor  of  10. 

5.  Light  Shift 

The  light  shift  was  investigated  for  both  cw  and  pulsed  operation  by 

using  two  different  Rb  lamps  (each  of  which  gives  a different  light  shift 

for  cw  operation.)  When  the  pulsed  optical  pumping  apparatus  is  operated 

-9 

cw,  a frequency  shift  (light  shift)  of  +6.5  x 10  occurs  when  the  lamp  is 
changed.  For  pulsed  operation,  the  frequency  shift  is  -2  x 10"10  under 
the  same  conditions.  It  is  unlikely  that  this  residual  frequency  shift  of 
-2  x 10  10  in  the  pulsed  mode  is  a light  shift  because  it  has  the  opposite 
sign,  but  even  if  it  were,  the  pulsed  mode  would  still  provide  more  than  a 
factor  of  30  reduction. 

This  frequency  shift  in  the  pulsed  mode  has  been  investigated  and 
found  to  be  proportional  to  light  intensity,  but  various  experiments  that 
we  have  carried  out  show  that  this  shift  does  not  behave  like  a true  light 
shift.  In  our  opinion,  this  shift  is  not  a true  light  shift,  and  hence- 
forth we  therefore  refer  to  it  as  the  "pseudo  light-shift  effect." 

The  pseudo  light-shift  effect  is  not  just  of  academic  interest.  It 

is  large  enough  that  a 1%  change  in  light  intensity  produces  a frequency 

-12 

shift  of  approximately  6 x 10  . It  is  therefore  of  great  importance  to 

understand  the  origin  of  this  effect,  and  to  reduce  (or  eliminate)  it. 

The  likelihood  of  accomplishing  this  is  considered  below  (see  conclusion 
7,  "Position  Shift  Effect".) 
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. Dependence  of  Frequency  on  Pulse  Parameters 


It  is  found  that  if  the  values  of  the  pulse  parameters  of  the 
pulsed  pumping  frequency  standard  are  changed,  the  output  frequency 
changes.  For  each  pulse  parameter  investigated  (T^,  T , T^).  it  is 
found  that  there  is  a range  of  values  where  the  frequency  changes  quite 
rapidly,  and  another  (separate)  range  of  values  where  the  frequency  chang- 
es much  more  slowly.  In  general,  it  has  been  found  that  it  is  possible  to 
choose  parameter  values  so  that  the  sensitivity  of  the  output  frequency  to 
small  changes  in  parameter  values  is  minimized.  In  general,  this  can  be 

done  without  having  to  sacrifice  signal.  This  approach  appears  adequate 

12 

for  stabilities  of  a few  parts  in  10  over  days  and  weeks,  but  if  stabili- 
13  14 

ties  of  parts  in  10  and  10  are  to  be  realized  over  the  same  periods  of 
time  (and  longer),  then  it  will  be  necessary  to  use  a more  fundamental  ap- 
proach. 

7.  Position-Shift  Effect 

The  "position-shift  effect"  occurs7  in  frequency  standards  of  the  type 

used  in  this  work.  It  is  due  to  inhomogeneous  broadening  in  the  physics 

package,  combined  with  spatially  nonuniform  optical  pumping  and  interroga- 
87 

tion  of  the  Rb  atoms.  In  general,  only  two  sources  of  inhomogeneous 
broadening  are  possible:  that  due  to  a nonuniform  static  magnetic  field 
(C-field  nonuniformity)  and  that  due  to  a true  light  shift. 

Based  on  the  NBS  findings,7  it  is  virtually  certain  that  the  pseudo 
light-shift  effect  is  due,  at  least  "n  part,  to  the  position-shift  effect. 
Indeed,  if  there  is  no  true  light  shift  in  the  pulsed  mode,  as  we  believe 
to  be  the  case,  then  the  position-shift  effect  could  be  completely  elimin- 
ated by  homogenizing  the  static  magnetic  field  (C-field);  this  would  re- 
move the  only  remaining  source  of  inhomogeneous  broadening.  Since  the 
static  field  in  our  physics  package  is  very  nonuniform,  its  homogeneity 
could  easily  be  improved,  and  the  position-shift  effect  significantly  re- 
duced. Therefore  it  might  be  possible,  by  using  this  simple  modification, 
to  greatly  reduce  the  pseudo  light-shift  effect.  This  modification  should 
also  reduce  the  sensitivity  of  the  output  frequency  to  changes  in  micro- 
wave  power.7 
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8.  Theoretical  Model  of  Device  Operation 


As  part  of  this  work,  a simple  theoretical  model  of  device  op- 
eration has  been  developed  that  provides  a simple,  conceptual  picture 
of  how  the  signal  at  resonance  depends  on  device  parameters,  especially 
the  pulse  parameters.  The  theoretical  predictions  of  this  model  are  found 
to  be  in  semiquantitative  agreement  with  the  experimental  results.  With 
further  work  it  should  be  possible  to  extend  this  theory  to  allow  predic- 
tions of  line  shapes  as  well.  The  importance  of  this  model  is  that  it 
provides  a conceptual  basis  for  understanding  the  physics  of  device  op- 
eration, which  is  important  for  a fundamental  approach  to  device  improve- 
ment. 

9.  Relaxation  Times 


A knowledge  of  relaxation  times  is  also  important  for  a fundamental 
understanding  of  device  operation.  One  of  the  results  of  this  work  has 
been  the  development  of  a new  method  for  the  measurement  of  the  relaxation 

time  T.  (time  constant  for  decay  of  the  population  difference.)  This  meth- 

1 33  34 

of  is  complementary  to  previous  methods.  * 

Measurements  of  Tj  for  our  experimental  conditions  show  that  the  dom- 
inant relaxation  mechanism  is  spin-exchange.  Our  best  set  of  measurements 
provides  evidence  for  the  presence  of  two  distinct  relaxation  times.  While 

this  might  be  due  to  a spurious  instrumental  effect,  it  could  also  be  due 

87  87 

to  two  differing  values  of  T.  for  our  apparatus:  one  for  Rb-  Rb  spin- 

* 87  85 

exchange  relaxation,  and  the  other  for  Rb-  Rb  spin-exchange  relaxation. 

If  this  result  is  real,  it  would  be  of  considerable  interest  because  the 
87  85 

Rb-  Rb  spin-exchange  relaxation  rate  has  not  been  previously  measured. 

10.  Line  Shapes 

The  rubidium  resonance  line  shapes  for  both  pulsed  and  cw  operation 
have  been  measured.  From  this,  the  following  conclusions  can  be  drawn. 

(1)  For  line  widths  greater  than  about  360  Hz,  the  line  width  is  propor- 
tional to  1/T^  where  T^  is  the  duration  of  the  microwave  pulse.  (2)  The 
"intrinsic"  line  width  for  pulsed  operation  in  our  apparatus  is  approxi- 
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mately  300  to  350  Hz  and  is  due  almost  entirely  to  spin-exchange  re- 
alxation.  (3)  For  cw  operation,  the  line  width  is  approximately  1100  Hz 
and  is  due  mostly  to  light  relaxation.  Therefore,  for  our  experimental 
conditions,  narrower  line  widths  can  be  obtained  when  the  apparatus  is 
operated  in  the  pulsed  mode.  (4)  The  line  shapes  for  pulsed  operation 
snow  small  subsidiary  maxima  in  the  wings  of  the  line.  This  type  of  line 
shape  has  also  been  observed  in  molecular  beam  experiments,  and  is  due  to 
the  pulsing  of  the  microwave  radiation. 

RECOMMENDATIONS 


1.  Long-Term  Stability  Measurements 

Phase  comparison  data  should  be  taken  for  the  pulsed  optical  pump- 
ing frequency  standard  over  a period  of  several  weeks,  or  more.  This 
will  allow  calculation  of  Oy(x)  for  x = > 24  hours.  Such  an  effort  is 
currently  in  progress  under  controlled  environmental  conditions  using  the 
house  frequency  standard  as  a reference.  These  measurements  will  most 

likely  be  limited  by  the  verifiable  (traceable  to  NBS)  stability  of  the 

-12  ,4 

house  standard  to  a precision  of  < 4 x 10  for  x = 24  hours  (8.6  x 10 

sec),  and  < 2 x 10  for  x = 1 week  (6  x 10  sec).  If  these  measurements 
yield  a Oy(x)  that  is  limited  by  the  stability  of  the  reference,  then  it 
will  be  necessaty  to  have  access  to  a more  stable  (for  x > 24  hours)  re- 
ference to  establish  the  true  value  of  Oy(x)  for  the  pulsed  optical  pump- 
ing frequency  standard.  Perhaps  0NR  could  arrange  for  a cesium  beam  stand- 
ard to  be  made  available  for  such  measurements.  Other  possibilities  may 
also  exist. 

2.  Pseudo  Light  Shift 

The  pseudo  light-shift  effect,  as  it  presently  exists,  is  expected  to 
limit  the  long-term  stability  (x  < 6 x 10  sec)  to  a few  parts  in  10  . 

For  this  reason,  it  is  important  t..jt  this  effect  be  reduced  (or  eliminat- 
ed). More  data  need  to  be  taken  to  establish  clearly  how  this  effect  de- 
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pends  on  device  parameters,  and  to  determine  its  physical  origin  (i.e., 
is  it  due  solely  to  the  position-shift  effect,  o*'  are  additional  physical 
mechanisms  involved?).  Also,  to  the  extent  that  it  is  due  to  the  position- 
shift  effect,  it  should  be  possible  to  reduce  it  significantly  by  improv- 
ing the  C-field  homogeneity  (which  is  easy  to  do).  Therefore,  this  ap- 
proach should  definitely  be  tried. 

3.  Position-Shift  Effect 


It  has  been  established7  that  the  position-shift  effect  is  responsi- 
ble for  frequency  changes  due  to  changes  in  light  intensity  (a  form  of 
pseudo  light  shift)  and  microwave  power.  Improvement  of  C-field  homogen- 
eity should  therefore  be  made  to  see  how  this  affects  the  pseudo  light 
shift. 

At  the  present  time,  no  investigation  of  the  effect  of  microwave  pow- 
er on  output  frequency  has  been  made.  This  should  be  examined,  at  least 
briefly.  If  there  is  a significant  dependence,  it  should  be  reduced  by 
homogenizing  the  C-field.  Another  approach  that  can  be  used  is  to  replace 
the  present  rf  circuit  board  by  one  that  gives  a more  amplitude  stable  out- 
put (such  a board  is  presently  available  to  us.) 

4.  Frequency  Shifts  due  to  Changes  in  Pulse  Parameters 

This  effect  will  affect  the  long-term  stability  at  the  level  of  parts 

13 

in  10  . Therefore,  understanding  this  effect  has  lower  priority  than 

understanding  the  pseudo  light-shift  effect.  Nevertheless,  an  understand- 
ing of  the  pulse  parameter  effect  is  desirable  because  it  is  likely  that 
there  are  physical  mechanisms  that  are  common  to  both  it  and  the  pseudo 
light-shift  effect.  Also,  it  is  unreasonable  in  practice  to  attempt  to 
separate  rigidly  these  two  effects.  All  of  these  various  effects  are  sure- 
ly interrelated,  and  complete  confidence  in  device  operation  will  require 
at  least  a rudimentary  understanding  of  their  origin. 

5.  Effect  of  Pulsing  on  Microwave  Spectrum 

38 

It  is  well  known  that  the  effect  of  pulsing  a sinusoidal  wave  is  to 
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generate  sidebands  that  are  placed  symmetrically  about  the  carrier  fre- 
quency. These  side  bands  are  equally  spaced,  with  the  spacing  between 
adjacent  frequency  components  being  equal  to  1/T  where  T = microwave 
pulse  duration.  The  amplitude  of  these  components  is  determined  by  the 
pulse  repetition  rate.  The  effect  of  these  additional  microwave  frequen- 
cies on  device  operation  should  be  examined  to  see,  for  example,  if  it  is 
possible  to  lock  the  atomic  resonance  to  one  of  these  components  rather 
than  to  the  carrier  frequency.  The  possibility  of  frequency  pulling  by 
sideband  asymmetry  should  be  studied  in  addition.  Also,  the  large  fre- 
quency shifts  that  occur  as  a function  of  T^,  when  T^  is  small,  (see  Fig- 
ure 12)  may  be  related  to  the  amplitude  and  location  of  these  sidebands 
in  some  way.  This  is  an  area  that  demands  investigation. 

6.  Temperature  Coefficient  of  Frequency 

At  the  present  time,  an  integrated  (resonance)  cell  is  being  fabrica- 
ted which  uses  a combination  of  buffer  gases  that  is  designed  to  reduce 
the  resonance  cell  temperature  coefficient  (TC).  This  should  demonstrate 
the  extent  to  which  the  resonance  cell  TC  can  be  compensated  using  physical 
mechanisms,  and  whether  this  approach  will  result  in  reduced  signal-to- 
noise  ratios. 

Stabilization  of  the  electronics  against  changes  due  to  variations  of 
ambient  temperature  is  basically  an  engineering  problem,  rather  than  a 
question  of  principle.  For  this  reason,  stabilization  should  be  carried 
out  only  to  the  extent  required  for  meaningful  long-term  stability  measure- 
ments. 

7.  Theoretical  Model  of  Device  Operation 

The  present  model  should  be  extended  only  to  the  extent  required  to 
aid  in  accomplishing  objectives  2 through  4,  above.  Emphasis  on  model  de- 
velopment should  be  placed  on  the  averaging  processes  that  occur  in  the  re- 
sonance cell,  and  on  understanding  the  observed  line  shapes.  For  under- 
standing device  operation  as  a frequency  standard,  the  theoretical  approach 
should  be  extended  to  include  FM  operation. 
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NOTE  ADDED  IN  PROOF 


- 12 

As  this  report  goes  to  press,  a value  of  Oy(x)  = 4 x 10  for  x = 

9 x 104  sec  (24  hours)  has  been  obtained  for  the  long-term  frequency  sta- 
bility of  the  pulsed  optical  pumping  frequency  standard  (see  below.) 

19 

a y(x)  is  the  square  root  of  the  Allan  pair  variance  for  four  samples  of 
fractional  frequency  (M  = 4),  each  having  an  averaging  time  x.  The  sta- 
tistical uncertainty  of  Oy  in  this  case  is  approximately  50  %. 

During  the  four  day  period  that  the  data  were  taken  the  maximum  vari- 
ation in  ambient  temperature  was  1.6  °C.  The  variation  in  the  average  am- 
bient temperature  during  successive  24  hour  periods  was  a few  tenths  of  a 
degree  centigrade  at  most.  Also,  the  unit  had  been  operating  for  more  than 

24  hours  when  the  data  taking  began.  Fractional  frequency  changes  due  to 

-12 

temperature  effects  should  therefore  amount  to  ~ 1 x 10  or  less. 

The  house  standard  was  used  as  the  frequency  reference  in  these  meas- 
urements and,  as  already  stated  in  this  report,  may  be  a limiting  factor  in 
the  value  of  Oy  that  has  been  obtained  above. 

It  should  be  noted  that  the  result  for  Oy  reported  here  serves  to  re- 
inforce the  conclusions  of  this  report. 
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APPENDIX  A 


THEORETICAL  MODEL  OF  DEVICE  OPERATION 


In  this  appendix  a theoretical  model  is  developed  that  can  be  used 
to  describe  device  operation.  At  the  present  time  this  model  is  re- 
stricted to  the  prediction  of  signal  strengths  as  a function  of  device 
parameters,  other  than  the  frequency  of  the  applied  microwave  radiation. 
That  is,  a theory  of  the  line  shape  is  not  included  at  the  present  time. 

This  model  begins  with  the  fundamental  quantum  mechanical  equations 
describing  the  optical  pumping  process,  and  subsequently  incorporates 
characteristics  of  the  present  device. 

Quantum  Mechanics  of  Pulsed  Optical  Pumping 

In  this  section  we  use  a simplified  version  of  the  theory  given  by 
34 

Vanier.  The  theory,  as  developed  in  this  section,  applies  to  an  opti- 

87 

cally  thin  layer  of  Rb  atoms  within  the  resonance  cell  when  the  applied 
microwave  field  is  resonant  with  the  rubidium  atoms.  Our  purpose  here  is 
to  try  to  obtain  closed-form  analytical  expressions  that  are  subject  to 

a simple  physical  interpretation  without  the  necessity  of  using  a compu- 
1 ter.  As  it  turns  out,  this  can  be  done  in  a meaningful  way  that  still 

allows  us  to  relate  the  end  product  to  the  actual  experimental  situation. 

We  begin  by  assuming  that  the  four  assumptions  stated  on  p.  129  of 
Vanier's  paper  are  satisfied.  While  this  is  not  strictly  true  in  prac- 
tice, the  deviations  from  actual  behavior  due  to  these  assumptions  should 

87 

, be  small.  Vanier's  theory  treats  the  Rb  atom  as  a nine-level  system: 

i eight  ground  state  sublevels  and  one  optically  excited  state.  Calcula- 

20 

tions  done  using  the  full  nine-levels  indicate  that  it  is  a good  ap- 
proximation to  consider  the  ground  state  as  composed  of  only  two  levels, 

1 and  2,  corresponding  respectively  to  (F=l,  Mp=0)  and  (F=2,  Mp=0).  This 
results  in  an  enormous  mathematical  simplification  (8x8  matrices  + 2x2 
matrices).  In  fact,  this  three-level  system  is  more  like  a two-level  sys- 
tem because  the  excited  state  does  not  appear  explicitly  in  the  formalism. 
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Following  Vanier's  notation,  the  equations  for  the  elements  of  the 
density  matrix  are,  for  a two-level  ground  state. 


p = -jrp-Vp-i-)  *'“■*)  (S1> 

fa,  = ^ T"*  p — ) (A2) 


" r +1^  fz  ~ «■  p12.  (A3) 

where  the  dot  indicates  the  partial  derivative  with  respect  to  time; 
where  the  notation  p = p^,  p2  = P22  has  been  used;  where  Tj  = 1/yj,  and 
= I/Y2  are  the  time  constants  for  the  decay  of  the  population  diff- 
erence and  the  oscillating  magnetization,  respectively;  where  r is  the 

87 

pumping  rate  due  to  the  light;  where  h>0/2tt  is  the  Rb  resonant  frequency 
for  the  F=2,  Mp=0»fcF=l,  Mp=0  hyperfine  transition;  where  io/2tt  is  the 
frequency  of  the  applied  microwave  radiation;  and  where  Aco^tt  is  the  true 
light  shift  that  is  present  only  when  the  light  is  on,  but  is  zero  when 
the  atoms  are  being  interrogated  by  the  microwave  field.  For  pulsed  pump- 
ing with  optical  detection  by  observing  the  absorption  of  the  pumping 
light,  the  important  quantity  is  p(t).  Since  the  pumping  light  and  the 
applied  microwave  radiation  are  not  present  simultaneously,  it  suffices 
to  solve  the  above  equations  for  three  special  cases.  This  can  easily  be 
done  using  standard  methods,  plus  the  results  given  in  Vanier's  paper. 

In  this  way  we  obtain 


p(rt  = Cpto)-j3  e"T|t 

(A4) 

pit}  = C + [pM-Cj  t 

(A5) 
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where 


pco  = i - Li  -pco>]  eT*  • 


• (cos  $t  + yl  sirvst) 

4C0  O 

(A6) 

S = V(*| if  -(Aj/2)1' 

(A7) 

c s U+r/^f1 

(A8) 

b s (r/z  + y() 

(A9) 

T = 

(A10) 

Ay  = X*-  y, 

no  light  or  microwave  radiation  present. 

(AH) 

Eq.  (A5)  is 

for  light  present,  but  no  microwave.  Eq.  (A6)  is  for  microwave  present, 
but  no  light.  In  each  case,  p(0)  is  the  initial  value  of  p for  the  con- 
ditions stated. 


AM  Signal 

The  quantity  that  is  observed  experimentally  is  the  AM  signal  at  re- 
sonance, denoted  by  S^.  How  this  signal  is  generated  is  discussed  in 
the  section  entitled,  "AM  Configuration".  Here,  we  relate  SftM  to  p(t). 

The  light  absorption  at  time  t is  proportional  to  p(t)  and  the  inten- 
sity of  the  pumping  light,  so  that 


absor  ption  OO  — IcF  jilt) 


(A12) 
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where  k is  a constant.  The  total  light  absorption  during  the  duration 
of  the  light  pulse  is  therefore  equal  to  the  integral  of  Eq.  (A12)  from 
0 to  T^.  The  average  light  absorption  over  one  cycle  of  pulsing  is  there- 
fore 

T* 

a = absorption(t)  = (T*  + 1 kT  dt  (Ai3) 

o 

The  observed  signal  is  the  difference  between  the  values  of  a with  the 
microwave  on  (6  t 0)  and  off  (B  = 0),  s^M  = a(s)  - a(0),  where  the 
lower  case  s indicates  that  this  is  the  contribution  to  S.^  from  an  op- 
tically-thin  layer  of  Rb  atoms  in  the  cell.  It  is  convenient  to  write 
this  as 


SAM  = A«AM/Oi  + Td)  <««> 

A aAM  = (i-e'b1*)  (ai5) 

where  p ^ is  the  initial  value  of  p at  the  beginning  of  the  light  pulse, 

a 

during  that  portion  of  the  10  Hz  cycle  when  the  microwave  radiation  is  not 
gated  off  (light  and  microwave  pulses  alternated).  Similarly,  p ^ is  the 

a 

same  quantity  when  the  microwave  radiation  is  gated  off  (light  pulses  only, 
and  no  microwave  pulses).  In  effect,  Aa  is  the  total  change  in  the  aver- 
age absorption  between  microwave  on  and  off  with  the  duty  cycle  factor  (T 

- 1 ^ 

+ Td)  removed.  Finally,  it  should  be  noted  that  Eq.  (A5)  was  used  in 

Eq.  (A13)  to  obtain  Eqs.  (A14)  and  (A15). 


75 


Steady-State  Behavior 


, Figure  A1  shows  (schematically ) the  time  behavior  of  p for  one 
cycle  of  pulsing.  In  the  treatment  given  here,  Ts  can  be  assumed  to  be 
zero  without  loss  of  generality. 

In  what  follows,  please  keep  in  mind:  that  the  population  of  the  F-l 
level  is  proportional  to  p(t);  that  the  value  of  p is  % when  the  system  is 
allowed  to  come  to  equilibrium  in  the  dark  with  no  microwave  present;  that 
the  effect  of  the  light  is  to  produce  a population  inversion  [p<h)  by 
pumping  atoms  from  the  F=1  level  to  the  F=2  level;  that  C is  the  value  of 
p when  the  light  has  been  present  for  a very  long  time  with  the  microwave 
not  present  (i.e.,  C = value  of  p for  cw  lamp  operation  with  microwave 
off).  Under  steady-state  pulsing  conditions  we  would  expect  behavior  of 
the  type  shown  in  Figure  Al:  At  the  beginning  of  the  light  pulse,  p has 
the  value  p , which  is  close  to  h (slight  population  inversion).  The  ef- 

a 

feet  of  the  light  is  to  pump  atoms  out  of  level  1 (F=l)  to  level  2 (F=2) 
so  that  the  population  of  level  1 decreases  exponentially  in  accordance 
with  Eq.  (A5)  (with  p(0)  replaced  by  p ) toward  the  asymptotic  (cw)  value 

a 

of  C.  At  point  b,  p = p^,  the  light  goes  off,  and  the  microwave  goes  on. 
The  net  effect  of  the  microwave  is  now  to  transfer  atoms  from  level2  to 
level  1.  The  population  of  level  1 therefore  increases  rapidly,  accord- 
ing to  Eq.  (A6).  At  point  c,  p = pc,  and  the  microwave  goes  off  leaving 
the  atoms  to  slowly  relax  undisturbed  in  the  dark.  The  population  of  level 
1 during  this  portion  of  the  pulsing  cycle  therefore  asymptotically  ap- 
proaches the  equilibrium  value  of  according  to  Eq.  (A4).  At  the  end  of 

the  pulsing  cycle,  p = p < The  condition  for  steady-state  operation  is 

a 

that  p = p7.  Self  consistency  also  requires  that 

a a 

-bit 

Pb  = = C + (pa  - C)  C (A16) 

. . -tTm 

Pc  “ ~ 2 ~ (l  ~ pb)  C 

• (cos  •+•  Ys  Sl w (A17) 
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po  = pm( V)  = X - (i  -pc)  fc" ^ (A18> 

where  we  have  used  Eqs.  (A5),  (A6)  and  (A4)  for  pj,  pjj,  Pjjj,  respec- 
tively. Eqs.  (A16),  (A17)  and  (A18)  constitute  three  equations  in  the 
three  unknowns  p , p.  , p in  terms  of  (supposedly)  known  physical  para- 

cl  D C 

meters.  For  our  purposes,  however,  it  suffices  to  solve  only  for  p.  and 

a 

to  use  the  result  in  Eq.  (A15).  In  fact,  if  Eq.  (A15)  is  to  be  used,  it 

(8) 

is  necessary  to  have  two  values  of  p . The  first  is  p ' ' which  has  al- 

3 (0)  a 

ready  been  discussed.  The  second  is  p ' ' which  can  be  obtained  by  sol/- 

a 

ing  for  pfl  in  Eqs.  (A16)  to  (A18)  with  8=0  (and  therefore,  s = i • Ay/ 2 , 

Tp  " 0.  and  Td  = Td  )* 

We  have  carried  out  this  solution  in  the  general  case  but  the  result 
is  very  complicated  and  not  very  illuminating.  For  this  reason,  we  will 
now  make  one  simplifying  assumption,  namely,  that  the  lamp  is  pulsed  on  for 
a long  enough  time  that  p approaches  the  value  for  cw  lamp  operation  (i.e., 
p^  = C).  If  this  assumption  is  not  made,  then  pfa  becomes  a complicated 
function  of  the  device  parameters  and  the  initial  condition  for  the  period 
following  the  end  of  the  light  pulse  depends  on  what  happens  during  each 
part  of  the  cycle. 

Mathematically,  the  condition  that  T^  -*  °°  is  equivalent  to  the  above 
assumption.  With  this  assumption,  Eq.  (A16)  becomes  p^  = C.  If  Eqs 
(A17),  (A18)  and  the  equation  p.  = C are  now  solved  simultaneously  for 

(o\  (01  D 

pfl,  then  both  and  p ' ' can  be  obtained  and  the  results  substituted 

into  Eq.  (A15).  If  we  denote  AaAM  as  T^  -*•  °°  by  AaAM,  then 


A~AMsUm  /AoAM)  ( 

Ao*m=  Jjp (j - C) eT|Ta  • 

• [l  - C ( cos  sTu  + —■  sin  sTu)] 


(A20) 
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Since  3 ~ 1600  sec'^  and  Ay  - 200  sec  * for  our  experimental  apparatus, 
it  is  an  excellent  approximation  to  write  s = 26  and  thus  A', /2s  = Ay/46 
~ 0.03.  This  allows  Eq.  (A20)  to  be  simplified. 


aoam  = *£■  (£-C)  e'¥,T<l  • 


[i  - e cos(2pT^l] 


(A21 ) 


Average  over  the  Resonance  Cell 


Eq.  (A21)  is  the  contribution  to  the  AM  signal  (aside  from  the  duty 
cycle  factor)  from  an  optical ly-thin  slice  of  the  resonance  cell  (see  Fig- 
ure A2).  It  now  remains  to  sum  the  contributions  from  all  such  slices. 
Before  doing  this,  we  note  that  the  pumping  rate  is  a function  of  z so 
that  b and  C are  also  functions  of  z according  to  Eqs.  (A9)  and  (A8).  We 
can  therefore  write 


I 


L 


L 

J AoAMU>  Az 


o 


(A22) 
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... 


1 ur*  A2 . Schematic  diagram  of  rubidium  resonance  cell  (not  to 

scale) . 


80 


AAam 


_L  c“*.Ta  . 

L C 


• J3fe)(l-e  ^ COS  [2.|bCX,>T^  ]]  dz. 


( A2  3 ) 


A/ 


Sam  " *AAM/(TA+Td) 


(A24) 


9CZ>S  ^ f^"C£z)J  (A25) 

Before  we  can  proceed  further,  we  need  to  know  the  functional  form  of  the 
pumping  rate  T[ z),  and  the  microwave  field  parameter  B(z). 

Since  the  resonance  cell  fills  essentially  the  entire  microwave  cav- 
ity in  our  apparatus,  the  microwave  field  in  the  resonance  cell  is  (to  a 
first  approximation)  the  same  as  in  the  cavity.  Since  the  cavity  operates 
in  the  TE^.,  mode,  the  z-component  of  microwave  magnetic  field  is  given 
by  Hz(z)  = Hq  sin  (ttz/L).  We  can  therefore  write 

!0  sin  (7VZ/L)  (A26) 

where  B0  is  an  average  value  of  3 over  the  disk  at  the  center  of  the  cell. 

The  integral  in  Eq.  (A23)  still  cannot  be  evaluated  because  r(z)  is 
unknown.  Even  if  r(z)  were  known,  the  resulting  expression  for  g(z) 
would  probably  be  sufficiently  complicated  that  the  integral  in  Eq.  (A23) 
could  not  be  evaluated  analytically.  For  this  reason  we  consider  two  sim- 
ple special  cases,  each  of  which  contains  some  element  of  the  truth. 
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Case  (1):  g(z)  = constant  = gQ 

This  case  corresponds  to  a uniform  optical  pumping  rate  over  the 
entire  resonance  cell.  The  integral  in  Eq.  (A23)  can  now  be  evaluated 
in  closed  form.  The  result  is 

AAAM  = g0  e T>  a fl  - £ x T|"  • 0o(apoTf>')]  (A27) 


where  Jq ( x ) is  the  zeroth  order  Bessel  function  with  argument  x. 

Eq.  (A27)  has  been  plotted  in  Figure  18  as  the  dash-dot  curve.  The 
constants  gQ  and  FQ  have  been  treated  as  adjustable  parameters  and  their 
values  chosen  to  give  the  best  fit. 

Case  (2):  g(z)  = step  function 

Figure  A3  shows  the  form  assumed  for  the  excitation  function  g(z),  in 
this  case.  In  the  figure,  g.  is  an  adjustable  constant.  For  this  choice 

i 7 

of  g(z),  the  rear  of  the  cell  is  preferentially  pumped  by  the  light. 

When  this  functional  form  of  g(z)  is  substituted  into  Eq.  (A23),  it  is 
possible  to  evaluate  analytically  the  integral.  The  result  is  an  infinite 
series  of  Bessel  functions  of  even  order.  In  this  case  it  is  necessary 
to  carry  out  numerical  calculations  through  J^x).  This  has  been  done 
and  the  results  are  presented  in  Figure  18  as  the  two  dashed  curves.  (The 
two  curves  are  for  different  values  of  the  relaxation  times.)  Additional 
discussion  of  the  results  obtained  here  can  be  found  in  the  main  body  of 
this  report  in  the  section  entitled,  "AM  Signal  as  a Function  of  Pulse 
Parameters. " 


